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ABSTRACT 


Broken rotor bars are a common cause of induction motor 
failures. In the past, the detection of broken rotor bars 
has primarily been limited to non-operating, and typically, 
disassembled machines. The ability to detect broken rotor 
bars while the machine is operating at normal speed and load 
is desirable. In support of the ongoing development of a 
failure analysis system for electrical machines, this thesis 
evaluates the method of using stator currents and voltages 
to detect the presence of broken rotor bars in squirrel-cage 
induction motors. The hypothesis of this method is that, 
given а sinusoidally applied voltage, the presence of 
certain harmonics in the stator currents could be used to 
detect the presence of broken rotor bars. i 


To support the evaluation, a system of first-order 
differential equations describing the electrical performance 
of a three-phase, squirrel-cage induction motor Was 
developed using stator phase currents and rotor loop 
currents as state variables. A FORTRAN simulation program 
was developed to solve the system of equations for a motor 
with or without a broken rotor bar. Using an "off-the-shelf" 
3-HP motor, numerical and physical experiments were 
conducted to test the failure detection hypothesis. Although 
the results of the numerical experiments indicated that the 
hypothesis was plausible, the experimental results showed 
that distinguishing between a manufacturing asymmetry and a 
broken rotor bar was impossible due to the existence of 
inter-bar currents. 


The existence of inter-bar currents in a squirrel-cage 
induction motor of the type tested in this thesis 
effectively "masks" the effects of a broken rotor bar. Thus, 
detection of broken rotor bars based upon a technique using 
only stator current measurements appears highly improbable. 
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CHAPTER 1 


INTRODUCTION 


1.1 Thesis Objectives 


In recent surveys 12,53) on the reliability of 
electric motors, three kinds of faults have been identified 
as constituting the majority of failures in induction 
motors. These are bearing-related (41%), stator-related 
(37%), and rotor-related (10%). Table 1-1 shows a summary 
of failures for these areas by components. The remaining 
failures (12%) are scattered among a variety of effects. As 
shown in Table 1-1, of the rotor-related failures, cage 
faults in the form of broken rotor bars or end rings are 
the cause of half of the machine failures. Cage faults 
occur due to design and manufacturing defects, misoperation 
and misapplication of the machine, lack of preventive 
maintenance, and aging/fatigue failure. Rotor cage defects 
result in machine failure due to increased frame 
vibrations, localized temperature increases on the rotor, 
and the "domino" effect in which one broken bar leads to 
another broken bar and so on. 

The purpose of this thesis is to present the results of 
research investigating the detection of broken rotor bars 
using stator current and voltage measurements. This 


research effort supports the ongoing development of a 
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failure analysis system for electric machines by the M.I.T. 
Laboratory for Electromagnetic and Electronic Systems 


[4,5]. 


Bearing-Related Stator-Related Rotor-Related 


(41%) (37%) (10%) 


Sleeve 16% Ground 
Bearings Insulation 


Anti-friction 8% Turn 
Bearings Insulation 


Seals 6% Bracing 


Thrust 5% Wedges 
Bearing 


Oil Leakage 3% Frame 
Other 3% Core 


Other 





. Table 1-1. Motor Failures by Components [1]. 


1.2 Condition Monitoring of Electric Machines 


The unexpected and sometimes catastrophic failure of 
electric machines can result in the reduction or total loss 
of production and operational safety, expensive repairs and 
extended downtime, and in most cases, large capital losses. 
In the case of military applications, these failures can 
result in the degradation of mission effectiveness and 


perhaps even the inability to perform a primary mission. 


11 


For many years, private industry and the military have used 
planned maintenance strategies to minimize electric machine 
failures. One major drawback of this strategy is that the 
need for corrective maintenance cannot be determined 
without removing the machine from service, disassembling 
it, and inspecting it. Without a method to externally 
determine the condition of an operating machine, a machine 
in perfect condition may be removed from service while a 
machine on the verge of failure maybe ignored. Obviously, a 
more efficient and cost-effective maintenance strategy 
would be to schedule maintenance and repairs based on a 
continuous assessment of the machine's condition while it 
is operating at normal speed and load. 

There are currently two methods being sed and/or 
developed to assess the condition of an operating electric 
machine [6]. The first method involves the analysis of 
vibration data and historical performance records. This 
method is commonly referred to as signature trend analysis. 
In signature trend analysis, sensor measurements (typically 
accelerations) are collected and processed through Fourier 
Transform at regular intervals. Each data collection is 
compared to previous data and known "good" baseline data in 
order to expose signature trends. Based on experience, 
these signature trends can be related to specific defects 
and failures. Thus, signature trend analysis is essentially 


a heuristic method used to assess the condition of a 
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machine. Although signature trend analysis is a viable and 
widely used monitoring system, the requirement for a 
database of historical performance and the experience 
required to relate a signature trend to a specific defect 
or failure are major diSadvantages. The second method is 
more theoretical by nature and seeks to identify the 
fundamental causes of failures and predict their effects. 
In this method, a detailed understanding of the machine 
response for various operating scenarios including normal 
operation as well as for various fault conditions 15 
required. Thus, for any given operating condition, the 
theoretical response of the machine 15 calculated and 
compared to the actual response of the machine. By 
comparing these responses an estimate or prediction of the 
machine's condition can be made. Although both of these 
methods typically require the use of a processor such as a 
er for continuous monitoring, each can be used 
both to prevent catastrophic failures by early detection 
and to aid in the preparation of routine maintenance 
schedules. 

A failure analysis system for electrical machines is 
currently being developed by the Laboratory of Electromag- 
netic and Electronic Systems (LEES) at the Massachusetts 
Institute of Technology. The development of this system is 
a major task of the Ship Service Power System Development 


research sponsored by the Department of the Navy’s Naval 


13 





Sea Systems Command [5]. The ultimate goal of the proposed 
failure analysis system is to provide a tool which can be 
used onboard naval ships to prevent, predict, and detect 
electric machine failures, and to suggest the corresponding 
maintenance. The system is intended for retrofit onto 
existing electrical machines or incorporation into future 
electrical machines and for implementation with inexpensive 
sensors and processors [4]. The following paragraphs, 
extracted from the research proposal [5], briefly describe 
the proposed failure analysis system. The system will use a 
combination of both methods described above to assess the 
condition of a machine. 

The underlying principle of the system is that 
failure prevention, prediction, and detection 
should be based on the estimation of states and 
parameters in relevant physical models of electri- 
cal machines. These models should include failure 
modes, mechanisms, or symptoms that are expressed 
in terms of the electrical machine states and 
parameters. The models, coupled with state and 
parameter estimation, then provide a means of 
directly and justifiably connecting measured data 
to impending or existing failures. Thus, the 
reliability of the failure analysis is enhanced. 

The models also provide a means by which measured 


data from a variety of sensors can be processed 
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together in a consistent manner. 

A (one-line) diagram of the proposed failure 
analysis system is shown in Figure 1-1. Measure- 
ments from sensors on the electrical machine are 
processed by a model reference state and parameter 
estimator. State and parameter estimates are then 


passed to a rule based evaluator which suggest the 


corresponding maintenance. Initially, the system 
will be developed using information readily 
available from terminal voltage and current 
sensors. The system can be expanded later to 


include information from other sensors such as 
thermocouples, accelerometers, acoustic sensors, 
and gas analyzers. 
The development of the failure analysis system 
is broken down into four subtasks. These are: 
1. Develop physical models which include 
failure modes, mechanisms, and symptoms 
expressed in terms of the model states and 
parameters. 
2. Develop estimators for the model states 
and parameters. 
3. Develop state and parameter evaluators 
that act on the estimated quantities so as to 
prevent, predict, or detect electrical 


machine failures and suggest corresponding 


15 





maintenance. 


4. Demonstrate and evaluate the results of 


theoretical work through numerical and 


physical experimentation. 


The results presented in this thesis directly support 


the physical model development of subtask (1) 


particular, this thesis is concerned with the 
and demonstration of models which include 


failures in induction motors. 


Electrical 
Machine 


Bensors 







Measurements 


State and 
Parameter 
Estimator 


Electrical 
Machine 
Model 


Estimates 





| Rule Based 


| te ы: State and 
Parameter 


Evaluator 


Suggested Maintenance 





above. In 
development 


rotor bar 


Evaluation and 


-- ----- -- Maintenance 
Rules 


Figure 1-1. Proposed failure analysis system [5]. 
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1.3 Detection of Rotor Bar Defects 


— — — 


In order to detect broken rotor bars, several methods 
have been employed. Visual inspection and bench test 
methods such as the "growler" and related probe techniques 
applied to disassembled motors have been used for many 
years [7]. The "single phase" test has also been used as a 
standard test for assembled but non-operating motors [7]. 
Recently, there have been a number of studies [7,8,9,10,11] 
to develop theories of the response of induction motors in 
the presence of broken rotor bars and/or end rings. 

A common result found in each of these studies is the 
existence of a lower sideband frequency component in the 
stator phase current when the motor is driven by a single 
harmonic stator voltage. This component of the stator 
current, which is at a frequency of (1-2s)f, where s is the 
rotor slip and f is the line frequency, is a result of the 
mime cai component of the backwards-rotating airgap 
field produced by the induced rotor currents. This field, 
which rotates backwards at slip speed with respect to the 
rotor, rotates forward with respect to the stator at 
(1-2s)Ns, where Ns is the synchronous speed. Thus, this 
field induces currents of frequency (1-2s)f in the stator 
windings. This component of the stator current causes 
torque pulsations and speed oscillations at twice the slip 


frequency. In addition to these effects, it has been shown 
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by Kliman et al. [7] and Penman et al. [11] that axially 
directed fluxes arise due to the asymmetry of the rotor 
magnetic circuit with a broken rotor bar. 

Based on these results, there have been a number of 
instruments developed to give an indication of a broken 
rotor bar while the motor is running at normal speed and 
load. These instruments detect broken rotor bars using 
measurements of one or more of the following parameters: 

- stator current 

- mechanical speed 

- frame vibration (i.e., acceleration) 

- air gap flux 

- axial leakage flux 
The most successful of these instruments has been an 
instrument based on detecting the twice slip frequency 
speed oscillations via a shaft position measurement [8]. 
However the sensitivity of the speed variation method is 
highly dependent on knowledge of the load inertia and 
torque, and results may be confused by other asymmetries. 

Recently, two new instruments have been used to 
successfully detect broken rotor bars in operating 
induction machines. Although both instruments are still in 
the "field test" stage of development, they appear quite 
promising. The first instrument, developed by Kliman et al. 
[7], uses stator current measurements to detect broken 


rotor bars. The second instrument, developed by Penman et 
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а1. [11], uses axial flux measurements to detect broken 
rotor bars as well as several other fault conditions such 
as phase imbalances and short-circuited turns on a stator 
winding. The following paragraphs summarize the operating 


theory and experimental results for these new instruments. 


1.3.1 Broken Bar Detector (BBD) 


The BBD being developed by Kliman et al. is based on 
existing theories for the performance of induction motors 
with broken bars. A broken rotor bar is modeled by 
superimposing a fault current, a current equal and 
Opposite to the normal current, on the rotor bar. The 
Magnetic field in the airgap caused by the fault current 
is always two-pole and rich in harmonics. This field 
anomaly rotates at the mechanical frequency of the rotor 
Since it is attached to the broken bar. The field can be 
resolved, using Fourier analysis, into an infinite series 
of counter-rotating, slip-frequency waves of smaller and 


smaller wavelength on the rotor. In addition to broken 


rotor bars, other fault conditions such as cage 
misalignment (rotor out of round), bearing misalignment 
(rotating eccentricity), and non-uniform magnetic orienta- 


tion of the rotor laminations create airgap field 
anomalies with a fundamental component on the same order 
as that of a broken bar. However, the higher-order 
harmonic components of the airgap field due to these other 


fault conditions are predicted to be much smaller than 
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those due to a broken bar. Figure 1-2 (from [7]) shows a 
comparison of the airgap flux spectra predicted for these 


various asymmetries and a broken rotor bar. 


0 1. Broken rotor bar 
2 Rotating eccentncity 
-10 3 Magnetic orientation in rotor laminations 


4 Rotating ovality (rotor out of round) 


- 80 


– 90 





- 100 


0123 45 6 7 8 9 10 11 12 13 14 15 16 
Air Gap Harmonic 


Figure 1-2. Airgap flux spectra for asymmetries [7]. 
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Based on this spectra, two conclusions can be drawn. 
First, with sensitivity sufficient to detect a broken bar, 
these other asymmetries may give rise to a false broken 
bar indication. Second, by examination of the higher order 
harmonic amplitudes, broken bars can be distinguished from 
asymmetries. The analytical expression for the frequencies 


present in the airgap flux is [7] 
tt (aaa! (1=1ў 
р 3 


where 


f,=line frequency 
k = harmonic index (1,2,3...) 
p= number of pole-pairs 


s=rotor slip 


For a typical motor, due to the design of the stator 
windings, only the odd, non-triplet harmonics of airgap 
flux couple with the stator windings. Thus, only those 
harmonic frequencies where k/p is 1,5,7,11,13, etc. appear 
in the stator currents. Table 1-2 (from [7]) shows these 
predictable stator current harmonics and relative 
amplitudes (as a percentage of the input current 


amplitude) due to an open rotor bar. 
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HARMONIC FREQUENCY AMPLITUDE 
(Hz) (%14.ш) | 

s | 
FUNDAMENTAL 


LSB 1 60(1-2s) 
USB 5 60(5-4s) 


USB 7 60(7-6s) 
LSB 7 60(7-8s) 


LSB = LOWER SIDEBAND USB = UPPER SIDEBAND | 


! 
| 
| 
| 
| 


1 


4 





Table 1-2. Current harmonics due to broken rotor bars [7]. 


Based on the above theory, Kliman et al. developed the 
computer-based instrument in Figure 1-3. The instrument 
performs two basic functions: signal processing and the 
implementation of а decision algorithm. The signal 
processing involves sampling the external leakage flux and 
current signals and transforming these time functions via 
the fast Fourier transform into frequency spectra. Since 
the harmonic frequencies of interest shown in Table 1-2 
are slip dependent, a precise measurement of the rotor 
speed is necessary. This is accomplished using an 
externally mounted coil which picks up a strong axial 


leakage flux from the rotor end ring at the slip frequency 
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(sf). Combining this with a measurement of the line 
frequency from the current signal, rotor speed is 
determined to within 0.2 rpm. The spectral windows of 
interest are computed and a narrow search around the 
predicted frequencies is performed. The amplitudes and 
frequencies of the broken bar signals are stored in 
computer memory, and if desired, can be printed out in a 
tabular form. The spectral windows can also be viewed on 
the computer screen. Figure 1-4 shows the experimental 
results (line current spectra in the vicinity of 60 Hz) 
presented in [7] for a test motor with varying degrees of 
fault. The increase in the amplitude of the lower sideband 


signal with increasing fault is quite evident. 


[Other inputs 
| (such as vibration) 







Current | Three 


transformer | phase. 
60 Hz 


Computer 
(inc. disk drives) 









Keyboard 


Figure 1-3. BBD one-line diagram [7]. 


23 









- 10 
- 20 
- 30 
- 40 
-50 
dB -60 
- 70 
- 80 
- 90 
- 100 
- 110 


- 120 
471 522 574 626 67.7 248 


Frequency (H2) 


No Broken Bars 
No offset 





- 10 
- 20 
- 30 
- 40 
-50 
ав -60 
-70 
-80 
-90 
– 100 
- 110 
- 120 





Опе Вгокеп Ваг 


47 1 52.2 574 626 677 72.9 
Frequency (Hz) 


Two Broken Bars 





471 52.2 574 62.6 72.9 
Frequency (Hz) 


Figure 1-4. BBD experimental results for s=0.01 Kar 
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Two decision algorithms may be implemented in the BBD. 
The first is a "trend" algorithm that recognizes the 
change from a healthy motor to a motor with a cage defect. 
Two sets of line-current frequency spectra are compared to 
detect any changes in the harmonic sidebands. If a 
Significant change (>20dB) occurs in the first harmonic 
lower sideband, the higher-harmonic sidebands are 
examined. If any of these harmonics have changed 
significantly, a broken bar is declared. Otherwise it is 
concluded than an asymmetry other than a broken bar is 
present. The second algorithm is a single-test diagnostic. 
Based on experimental data gathered from 23 power plant 
motors, "good" motors exhibit a first harmonic lower 
sideband of -60dB or less relative to the line frequency 
component. Thus if the first harmonic lower sideband is 
less than -60dB there is probably no fault. If the first 
harmonic sideband is greater than -50dB there ‘is probably 
a broken rotor bar. These thresholds will be updated, if 


necessary, based on the results of further field tests. 


1.3.2 Leakage Flux Detector 


The leakage flux monitoring system developed by Penman 
et al. [11] is based on theories describing the harmonic 
content of the axial flux produced by an electric machine 
for various fault conditions. Ideally, with no faults, an 
electric machine will not produce a net axial flux. 


However, due to small asymmetries in both the material and 
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geometry of a machine, a small, axial leakage field is 
produced and can be detected by an externally mounted 
coil. Under the assumption that fault conditions represent 
a gross change in the electric and/or magnetic circuit 
behavior of a machine, the faults can be identified by the 
effect these changes produce in the axial leakage field. 
In order to use this technique, a master table of all 
possible axial flux harmonic components for a given 
machine type must be generated. Once this is done, certain 
groups of these harmonics can be identified with various 
fault conditions. 

For an induction motor, the harmonic components of the 
airgap flux produced by a balanced 3-phase stator winding 


is given by 


B= B,cos(wt-p6)+ Bscos(wt+5p0) 


+ B,cos(wt-7p0)+B,,cos(wt+11p0)... (1-2) 


where 


p=number of pole pairs 


9= angle from a reference point on the stator 
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То determine the harmonic components of the induced 
rotor currents, a transformation to the rotor reference 
frame is necessary. The transformation is made using the 


relation 


BOIS C. 


9 1-3 
p ( ) 
where 
w(l-s) 


= Totor angular velocity 


Ф = angle between the stator and rotor at t=0 


Thus, in the rotor frame of reference, the airgap flux 


density is given by 


B= B,cos(swt-p4)+ Bscos((6-5s)wt+5pp) 


> B,cos((7s-6)wt-7p0)+B,,cos((12-11s)wt+11pé)... (1-4) 


The harmonic frequencies of the currents induced in the 
rotor bars by this field are slip dependent. The small 
axial leakage flux, which is produced by the eae 
currents and minor asymmetries in the machine, will 
contain these same frequency components. If a fault 


condition exists, the expressions for the airgap flux 
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density given above аге no longer valid. For example, a 
phase asymmetry on the stator winding will result in the 


airgap flux density 


B- ) B,[cos(wWt-np0)+cos(wWt+np0)] (1-5) 


n odd 


which, in the rotor frame of reference, becomes 


| 
pee. B,cos[(l=n(1-s)wtenpé] (1-6) 


n odd 


Thus, additional frequency components are present in the 
induced rotor currents and axial leakage flux. 

Using the same methodology, a table of axial flux 
frequencies (or harmonics) can be developed for various 
fault : conditions. Table 1-3 (from [11]) shows the 
predicted axial flux harmonics for a 4-pole induction 
motor at a slip of 0.02. In the table, stator and rotor 
faults are grouped as either symmetrical (related to a 
phase) or unsymmetrical (related to a pole). The table 
includes effects due to the fundamental (f1) and third 
harmonic (f3) components of the supply current. The 
numbers in each column represent the component of the 


airgap flux wave that produces the associated axial flux 
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harmonic (i.e., 1 is the forward traveling first harmonic 
and -3 is the backward traveling third harmonic 
component). 

By monitoring the spectral components of the axial 
leakage flux, a fault group can be identified by comparing 
the spectra to the master table. This is the philosophy 
incorporated into the leakage flux monitoring system. A 
one-line diagram of this system is shown in Figure 1-5. 
The transducer (a printed-circuit split coil) is 
externally mounted on a machine. The transducer signal 
(induced voltage by the leakage flux) is processed and 
transferred to the diagnostic unit. The diagnostic unit 
analyzes the spectra and identifies any fault conditions 
present using the master spectrum table. As stated 
previously, this system has successfully detected various 
faults imposed on a 4KW squirrel-cage induction motor. 
Figures 1-6, 1-7, 1-8 are samples of the experimental 
results presented in [11] for various fault conditions. 
Figure 1-6 is the axial flux spectra for a healthy motor. 
Figure 1-7 shows the large amplitude increase in the 
fundamental (50 Hz) component for a short-circuited stator 
winding. For a broken rotor bar, increased amplitudes in 
the fundamental as well as third (150 Hz) and fifth (250 


Hz) components are evident in Figure 1-8. 
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Number Frequency Frequency Stator harmonics Rotor harmonics 





factor ато ОХО Sn HN 
Hz Phase Pole Phase Pole 
На A tr E A Y 
бб 4 6 h G k # 
———————— —— CNN 
1 s 10 1 2 
2 3s 30 3 6 
3 (51 21582 215 7 =1 
4 (i351 23082 235 3 il 
5 (1+92 255 1 
6 ПИ 2275 5 1 
7 46-1 460 -1 =@ 
8 2-4 480 EX =2 
9 1 500 1 2 
10 25-1 520 1 2 
11 (95 – 3) 2 705 9 "5 
12 (5s - 312 725 5 ==) 
13 13-512 745 -1 3 
14 (3-352 765 3 3 
15 552 950 5 10 
16 ego 970 3 6 
V7 2=% 990 =] -2 
18 24 545 1010 1 2 
119) (115-/5/2 1195 11 -5 
20 ss 5 2 2 ES 7 -5 
21 153507517315 -3 5 
22 (5158512 1255 1 5 
23 65-3 1440 -3 -6 
24 45-3 1460 -3 -6 
25 SES 1480 3 & 
26 2) 1500 3 6 
27 113-712 1685 13 ез 7) 
28 СОБЕРУ 1705 9 =7 
29 17-5892 1725 == 7 
30 ES 1745 = {] 7 
31 75-4 1930 7 14 
32 55-4 1950 5 10 
33 4 - 3s 1970 213 -6 < 
34 А5 1990 =$ == 
55 (11655 29120082 1/25 15 -9 
36 Has 9 ngu 11 -9 
e I8—7s2 2215 = у 9 
38 19- 3612 2235 zi 9 
39 85-5 2420 EE le) 
40 65 -5 2440 = - 10 
41 5-45 2460 š 5 10 
42 2-25 2480 5 10 
43 1175-1112 2665 17 -11 
54 (138-1112 2685 13 БАП 
45 111-9s)2 2705 -9 11 
46 5 572 0627245 5 11 
47 9s- 6 2910 9 18 
48 7s-6 293.0 7 14 
49 6 - 5s 295.0 esc -10 
50 6 - 3s 2970 ==) -6 


Table 1-3. Axial flux harmonics [11]. 
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Figure 1-5. Leakage flux monitoring system [11]. 
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Figure 1-6. Healthy motor [11]. 
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Figure 1-7. Short-circuited stator winding turn [11]. 
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Figure 1-8. Broken rotor bar [11]. 
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1.4 Significance of Thesis 


From the discussions presented in the previous 
sections, one can see that broken rotor bars can be (or 
have been) detected using stator current measurements. 
Thus, the results presented in this thesis confirm existing 
theories and experimental results [7,8,9,10,11]). However, 
several contributions are made to the field of electri- 
cal-machine fault analysis. 

First, a more general approach is taken to develop the 
system of equations describing the electrical performance 
of an induction motor. If desired, the system of equations 
and computer simulation can be used to analyze other fault 
conditions as well as broken rotor bars (subject to the 
limitations resulting from the assumptions discussed in 
Chapter 2). Although the system of equations and computer 
simulation developed consider only fundamental space 
ен, both can be easily modified to include 
additional space harmonics. In order to analyze other fault 
conditions and/or include additional space harmonics, the 
expressions for the appropriate matrix elements must be 
modified. This flexibility is not possible using the method 
presented by Williamson and Smith [10], where the equations 
were derived for a specific fault condition using specific 
stator current harmonics. In order to analyze other fault 
conditions or to include additional space harmonics, an 


entire new set of equations must be developed. 
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Second, a relationship between the standard sin- 
gle-phase, equivalent-circuit model of an induction motor 
and the system of equations has been derived. With this 
relationship, the electrical parameters, for example rotor 
bar resistance, rotor bar leakage inductance, and stator 
phase leakage inductance needed to solve the system of 
equations can be easily calculated from the single-phase 
circuit model values. Instead of disassembling the machine 
and measuring these parameters or requesting the design 
data from the manufacturer, the standard no-load, 
locked-rotor, and DC tests can be used to determine the 
required parameters. 

Finally, the experimental results provide another "data 
point” or quantitative comparison between motors with and 
without broken rotor bars. Combination of this data with 
the previous data collected and presented by others may aid 
in determining threshold values for fault analysis system 
decision algorithms. The experimental results presented 
show that the existence of inter-bar currents (currents 
which flow through the rotor iron) in squirrel-cage 
induction motors mask the existence of a broken rotor bar 
and thus, severely limit the ability to detect a broken 
rotor bar using stator current measurements. This finding 
confirms the analysis of inter-bar currents presented by 


Kerszenbaum and Landy [9]. 
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1.5 Outline of Thesis 


The material presented in the following chapters is 
organized in the same fashion in which the research was 
conducted. Chapter 2 starts out with the development of the 
system of equations describing the electrical performance 
of an induction motor. The assumptions used to develop the 
equations are discussed in detail. In addition, the 
relationship between the single-phase, equivalent-circuit 
model and the system of equations is presented. The 
solution to the system of equations is addressed in Chapter 
3. The numerical technique and computer simulation program 
used to solve the equations are described. Simulation 
results (using the parameters of the experimental test 
motor) with and without a broken rotor bar are presented. 
In Chapter 4, the experimental results using an 
"off-the-shelf" 3-HP motor with and without a broken rotor 
bar a up csented, analyzed, and compared to the simulation 
results. Finally, Chapter 5 summarizes the results of the 
research, discusses the limitations of the analysis, and 


provides some recommendations for follow-on research. 
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CHAPTER 2 


THEORETICAL DEVELOPMENT OF SYSTEM 


2.1 Discussion of Approach 


In order to meet the primary objective of this research 
which is to detect broken rotor bars using stator current 
and voltage measurements, the system of equations 
describing the electrical performance of an induction motor 
is developed using stator phase currents and rotor loop 
Currents as state variables. Each stator phase and rotor 
loop is described in terms of resistances and inductances. 
The forcing functions for the system are the stator phase 
voltages. The rotor is described in terms of loops in order 
to facilitate the determination of stator and rotor 
coupling. 

Using classical field theory and a coupled-circuit 
viewpoint in which stator phases and rotor loops are 
regarded as circuit elements whose inductances depend on 
the angular position of the rotor, flux linkages are 
expressed in terms of stator phase and rotor loop currents 
and inductances. The stator phase voltages are related to 
these flux linkages using Faraday’s law. The result is a 
system of first-order differential equations cae the 
electrical performance of an induction motor. The form of 


this system of equations is 
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VOJ FAORI O) (2-1) 


where 


(v(t))- voltage vector 

[A(t)]= flux linkage matrix = [1(0(1))] [i(t)] 
[R]= resistance matrix 

[i(t)]= current vector 

[L(9(t))]= inductance matrix 


6(t)= rotor position 


This approach is quite general and can be used to 
describe the electrical performance of any induction motor, 
healthy or failed. In this analysis, a transformation 
matrix will not be used to eliminate the rotor-position-de- 
pendent elements of the inductance matrix. The use of a 
transformation matrix requires the transformation of rotor 
currents as well. Since a broken rotor bar will be 
simulated by setting the rotor bar current to zero, the use 


of & transformation matrix is undesirable. 


2.1.1 Assumptions 


Several assumptions are made in order to reduce the 
complexity of the analysis. Although each assumption 
reduces the generality of the analysis, the goal of 


understanding the effect of a broken rotor bar on stator 
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currents can still be achieved. It should be noted that 
each assumption can be relaxed and the system of equations 
can be developed using the method presented. The following 
paragraphs summarize the assumptions used and their 
consequences. 
1.) The stator windings are balanced and sinusoidally 
distributed. Thus, the airgap field produced by stator 
phase currents consists of a fundamental space 
component only. In addition, only the fundamental 
Space component of the airgap field produced by rotor 
currents will couple with the stator windings. 
2.) With the exception of a fault, the rotor cage is 
symmetrical. This dictates that the spacing between 
adjacent rotor bars, as well as the rotor radius, is 
constant over the surface of the rotor. In addition, 
rotor skew is neglected. 
3.) The rotor cage end rings are perfectly conducting 
disks. Thus, no axial flux is produced and the sum of 
the rotor loop currents is zero. Thus, end effects are 
neglected. 
4.) The rotor bars are insulated from the rotor iron. 
Current flows only in the rotor bars and thus, 
inter-bar currents which flow through the rotor iron 


laminations are neglected. 
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5.) The airgap dimension is small compared to the mean 
rotor radius. Thus, the airgap magnetic field is 
assumed to be only radially directed and to be only a 
function of azimuthal angle. 

6.) Only the fundamental space component of the airgap 
field produced by rotor currents will be considered. 
Higher order space harmonics will be neglected since 
the stator windings couple with the fundamental 
component only. 

7.) The mechanical speed of the rotor is constant. The 
system of equations describe the electrical perfor- 
mance of the motor and the mechanical dynamics are not 
included. Rotor speed will be considered an input to 
the system. 

8.) The magnetic backing material on the rotor and 
stator has infinite permeability. Thus, the magnetic 
field intensity is confined to the airgap region. In 


addition, saturation effects are neglected. 


2.2 Development of Equations for an N-bar Rotor 


Before proceeding with the development of the system of 
equations some background information is necessary. The 
following paragraphs describe the nomenclature used, the 
geometry and coordinate system, and the definition of a 


rotor loop. 
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Wherever possible, an attempt has been made to use 
"standard" symbols for parameters throughout this thesis. 
The following list of symbols is provided to avoid 
confusion which may arise due to the numerous variations of 
“standard” symbols. The appropriate MKS units are included 


in brackets following the definition of each symbol. 


List of Symbols 





г,0,2 = cylindrical coordinates for stator reference 
frame 

т',.0',2' = cylindrical coordinates for rotor reference 
frame 

B = magnetic flux density [Wb/m?] 

B, = radial component of flux density [Wb/m?] 

d - active length of machine [m] 

E = electric field intensity [V/m] 


= airgap length [m] 


H : 2 magnetic field intensity [A/m] 

H, = radial component of magnetic field intensity 
[A/m] 

[1] = electrical current vector [A] 


| = electrical current [A] 

= current in stator phase a,b, or c [A] 
= current in rotor loop n [A] 

= current in rotor bar n [A] 

JT = current density [A/m?] 

K = surface current density [A/m] 


axial component of surface current density [A/m] 


ч 
" 
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[L] 


an - Lia 


n... 


[A] 


А 


и, 


inductance matrix [H] 

leakage inductance of rotor bar n [H] 

leakage inductance of a stator phase winding [H] 
mutual inductance between windings m and n (H] 


mutual inductance coefficient for a rotor loop 
and stator phase winding [H] 


total number of turns of a stator phase winding 
[turns] 


stator winding azimuthal turns density [turns 
per unit azimuthal length] 


number of rotor bars 

number of machine pole-pairs 
resistance matrix [Ohms] 

effective resistance matrix [Ohms] 
mean rotor radius [m] 

resistance of rotor bar n [Ohms] 
resistance of a stator phase winding [Ohms] 
rotor slip 

time [sec] 

voltage vector [V] 

stator phase a,b, or c voltage [V] 
electrical angle [radians] 


angular displacement between stator and rotor 
reference frames at t=0 [radians] 


flux linkage matrix [Wb] 
stator phase winding a,b, or c flux linkage [Wb] 
rotor loop n flux linkage [Wb] 


permeability of free space [H/m] 
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w = stator electrical frequency [rad/sec] 
о, = rotor angular velocity [rad/sec] 
w, = rotor electrical frequency [rad/sec] 


The standard definitions for electrical angle, rotor 
slip, and rotor electrical frequency are shown below. 
Unless otherwise specified all angles will be given in 


mechanical degrees. 


electrical angle: д,= ре (2-2) 
rotor slip: wW- pw, 
g= —T (2-3) 
о 
rotor frequency: о), = во - ро) "sw (2-4) 


The coordinate system chosen to describe the stator and 
rotor is a cylindrical system. The stator is described 
using a fixed coordinate system г, 0.2). The rotor 15 
described using a moving coordinate system (r’,@’,z’) which 
travels at an angular velocity of w, with respect to the 
fixed coordinate system. Defining 0, ав the angular 
separation between fixed and moving systems at t=0, a point 
on the rotor in terms of the stator reference frame is 


described by 
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2-2 (2-5) 
DRDS (56 


STATOR 





ROTOR 


Figure 2-1. Coordinate system. 


As stated previously, the rotor is 2222: іп terms 
of loops. Figure 2-2 shows a section of the rotor cage. A 
rotor loop can be thought of as a current path which 
includes two adjacent rotor bars and the portions of the 
end rings connecting the bars. Thus for an Na-bar rotor, 
there are Nas rotor loops (Na adjacent single-turn coils). 
Since the rotor is symmetric, the separation between 


adjacent rotor bars is 2nz/N,, radians. 
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BAR 


END RING — Loto BARS 


Figure 2-2. Rotor loop description. 


Rotor loop current is defined as the current flowing in 
the single-turn coil (see Figure 2-3). The actual rotor bar 
current can be expressed in terms of rotor loop currents. 


For example, the current in rotor bar n is given by 


ein ins (250) 


Each rotor loop current produces an airgap flux which 
couples with each stator phase winding and every other 
rotor loop. In addition to this coupling, every rotor loop 
couples with the two adjacent rotor loops due to the 
resistance and leakage inductance of the common rotor bars 


as shown in Figure 2-3. 
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Figure 2-3. Coupling between adjacent rotor loops. 


In order to derive the system of equations, two steps 
are required. First, the flux linkage of each winding 
(stator, phases and rotor loops) must be expressed in terms 
of stator phase and rotor loop currents and inductances, 
Second, a relationship between the flux linkages and 
forcing functions (stator phase voltages) is required. The 
flux linkage of each winding will be determined using the 
principle of superposition. That is, the total flux linked 
by a winding is simply the sum of the flux linkages due: to 
each current flowing in the system. The flux linkages can 
be related to terminal voltages using Faraday's law. In 


mathematical terms for & system with m windings we have 
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NT is (2-7) 


d 
Moos B, (2-8) 


where 
A,=flux linkage of winding n 
Lam = mutual inductance between windings n and m 
ја = current in winding m 
Va= terminal voltage of winding n 
R, = resistance of winding п 


i, "current in winding n 


For an induction motor with Nrs, rotor bars, there are Nas 
+3 windings (Nx rotor loops plus 3 stator phase windings). 
Since there is no net axial flux, the sum of the rotor loop 
currents is zero and the system of equations can be reduced 
from Net+3 to Nr+2. The following sections detail the 
calculation of inductances for each winding. The 
calculations are broken down into three sections; stator 
relations (involving stator terms only), rotor relations 
(involving rotor terms only), and stator-rotor relations 


(coupling between the rotor and stator). 


2.2.1 Stator Relations 


Each stator phase winding is assumed to be represented 


by a sinusoidally distributed surface winding, separated 
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by 120 electrical degrees from the adjacent phases. For a 


p pole-pair machine, the stator winding azimuthal turns 


density is given by 


5 








И oq совре 
N, 2n 
n SE сөз рө- 21) (2590 
8 


N, 2n 
n Pos рө + = 
9 


Open-circuiting phases b and c and exciting phase a 
with a current I, a surface current is produced on the 


Stator according to 


$ 


K N 
ER 





cosp® (2-10) 


Choosing a contour that crosses the airgap at 0=% and 
0-6-n/p, the radial magnetic flux density can be 


determined using Ampere’s law; see Figure 2-4. Thus, 


f mats f [1-826 tuo M 


$7n/p 


Н,(ө)а-Н,(Ф*л/р)а= |, K,R d0=- 
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ѕіпрф (202) 


47 





STATOR P 
| nm ЕЗ 
с^ | 9 
| 
| | АІРСАР 
ROTOR 


Figure 2-4. Integration contour on stator. 


Using the symmetry relation H,(¢)=-H,(¢+a/p), 





Н,(ф)= – 


апр . 2M 
ар ( ) 


SN. р 
БИЕ o np (2-14) 


Since é is an arbitrary angle, the radial flux density 


in terms of 0 is 





N 
B,(0)= a sinp® Gas) 
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The flux linked by a full-pitch, single-turn coil can 
be calculated by integrating over the surface enclosed by 
the coil. Note that the surface enclosed by a coil for all 
flux linkage calculations is obtained by traversing the 
coil in the direction of positive current flow, i.e., 
current flows in the z direction at 6 апа іп the -z 
direction at é*mz/p. For the coordinate system being used, 
the normal component of this surface is in the -r 


direction. In this case, 


d 6-л/р 
Ф(%) -f п B-nda 
9 4 


МАРИ МЕГ. 
-af sinp9R,d9 (2= 716) 
, 2ар 





о 
> === ос 


Again, since é is an arbitrary angle, the flux linked by a 


coil extending from angle 0 to angle 0*n/p is 


ави 
e#(0)= ° ° T cospg (2—17) 
gp 


The total flux linked by the stator phase windings can 
be calculated by integrating the number of turns times 
this flux linkage over the winding surface. This results 


in 
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n/p 
he -pf n,(0)#(0)R,d80 


1, NidR,I 


4gp? (2-18) 
Sa/3p 
À, -pf n,(8)6(0)R,d0O 
24/3р 
авл 
a A (2-19) 
7a/3p 
de -pf п,(0)Ф(Ө)К, Ө 
4a/3p 
-u,n NidR,l 
as (2-20) 


8gp* 


From these flux linkages, the inductances of the stator 


phase windings due to a current in the phase a winding are 


TON @ 


ss > T 4gp? M [ (2 =2 1) 
А, -nu,tN*R,d 

Шығын” Баг (2-22) 
A. -u.nN?R.d 

M оғ” (2-23) 


89р? 


where Lẹ is the leakage inductance of phase a representing 
additional phase a flux linkages not accounted for by the 


space fundamental component of the airgap flux. 
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Defining the inductance due to the space fundamental 
component of the airgap flux produced by a stator phase 
current as the stator phase airgap inductance (L;), the 


following relations are obtained 


ПЕ То m E ни. (2-24) 
Ls 
ЕЕЕ КИЕ (2-25) 
мћеге 


ui aNiR,d 


(2-26) 
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2.2.2 Rotor Relations 


Applying a current I to rotor loop 1 and open-circuit- 
ing all other rotor loops, the radial magnetic flux 
density can be determined using Ampere's law and Gauss's 
law. From Ampere's law, it can be shown that the radial 
flux density is a constant (B) in the region 0<9’<?n/Na 
and that the radial flux density is & constant (B;) in the 
region 2n/N.,<9 <2n (see Figure 2-5). Using Gauss’s law and 
noting that there is no axial flux, a relation between 


these constants can be calculated. In this case 
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4 2n/ Nay 2a 
|| i B,R, do’ [ 5, R, a0 az-o (2:382) 
2 


АГ Ма 


puc oL JB; (2-28) 


Choosing a contour which crosses the airgap and 
encloses rotor bar 1 (see Figure 2-5), the radial flux 
density can be determined from Ampere’s law. Following 


this analysis, 


STATOR 





e 
0 E 
ne 
Figure 2-5. Integration contour on rotor. 
НИН (uc 
g 


Using B- 4,H and equations 2-28 and 2-29, B, can be found as 
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_- Бо (Nga 1) 


B, 
IN» 


Since each rotor loop is a single-turn coil, 


linkages can easily be determined. In general, 


4 2ла/Н,, 
ха" | f B,(0")(-R,)d0 dz 
0 2 


2(a-1)/Nas 


for rotor loop 1, this yields 


21 IN: - 1) 
А GN ir 


and for all other rotor loops this yields 


-2nu,dR,l 
To q N, 


(2-30) 
the flux 
(2-31) 
(2232) 
(2-12 


From these flux linkage values the rotor loop 


inductances are determined to be 


Хи 2ли„К, 4 (Му - 1) 
I gNës 


Es a 


ubi * Linz 


Ань -2ng,R,d 
= ] аМ% 





Lil 


Lin 
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(2-34) 


(2-35) 





А > ои Ка 
I ghz ` Lrb2 (2236) 





En = 


where Lin. and Lins are leakage inductances of the rotor bars 
accounting for flux linkages not accounted for by the 
airgap flux distribution. 


For the other rotor loops (n=3,4,5....Nra-1) 


L _ Am 2лиоК,д (2 om 
rarl I аһ:, ) 


From equations 2-34, 2-35, and 2-36 one can see how 
the leakage inductance of each rotor bar is associated 
with the two adjacent rotor loops. Thus, for each rotor 
loop, the leakage inductances of the two rotor bars which 
make up the rotor loop are included in the "self 
inductance" term. For the two adjacent rotor loops, the 
leakage inductance of the common rotor bar is included in 
the "mutual inductance” term. 

Defining the inductance due to the airgap flux 
produced by a rotor loop current as the rotor loop airgap 


inductance (L,), the relations 
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Lr 


а = Eres m»*n,nzxl 


are obtained where 


2ap, R d/Npg-1) 
М 


R 


2.2.3 Stator-rotor Mutual Relations 


From section 2.2.1, the radial flux density 


airgap.due to a current I in stator phase a is, 


- u, NI 
ae anne 


ODE 


(2-99) 
(239) 
(2-40) 
(2-41) 
(2242) 
in the 
(2-43) 


To determine the flux linked by a rotor loop, this 


flux density must be transformed to the rotor frame of 


reference. In the rotor reference frame the radial flux 


density is 
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кешене. ple + + ) 2- 44 
т sin Ett 0 

2gp E ( І 
The flux linked by a rotor loop is 


4 2nn/Nap 
| || B,(0°)(-R,)d0‘dz (2-45) 


л(в-1)/8,, 


да 1 ` ЕЗ | [Gre 
= sin| — |Jsinpn| —— * o, t* 6, 
gp” N Nap 


RB 


The mutual inductance between phase a and rotor loop n is 


Ara 4, N,R,d 


2n-1 | 
Pn = sig Jeino[ C327. 8, vo.) (2-46) 
l gp RB Мав 


© 


2. 


Similarly, for currents applied to stator phase b and c, 


the results are 


N.R, d 2n-1)x 2n 
Lu eius (P ais (o( G0 DT. i, ceo, - 27] Ben) 
gp 


N¿R,d -] n 
пат Ме sinh ZB Jein (p FEDE, wto.) 7) (2-48) 
gp Маз 3 


Using the relation L,,-l, and defining M as the 
constant coefficient of the mutual inductance between a 


stator phase and rotor loop, the relations 
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are obtained where 





2.2.4 Overall System of Equations 


(2-49) 
(2-50) 
(2-51) 
(2-52) 


Using the relations derived in the previous sections, 


the flux linkages for the system can 


matrix .form 


[А] = [1 ][1] 
where 
A, 
À, 
А. 
ГА > А, 
TIT. 


57 


be expressed in 


(2-53) 


(2-54) 





La La La, m 240 Гоги ре 
i: Ly i Lon 216% Lyn RB 
[L]= 
Ies IT I INS ноти 
і, 
i, 
i, 
[ШИЕ ін 


Tey RB 


The flux linkages are related to 
by Faraday’s law. For example, 


stator phase a, the following is obtained 


ЕШ” 
dt 





УН = 


Similarly, for rotor loop n (whose net voltage 


because it is shorted), Faraday’s law yields 


ЛАС 


(В rho * КОЕ Ra Ius Ка кеј = dt 
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(2-55) 


(2-56) 


the stator phase voltages 


applying Faraday’s law to 


(257) 
is zero 
(2-58) 





Defining a resistance matrix, [R], and combining the above 


results, the system of equations can be expressed as 


i all 
[v}= (Li) (qm 1 Í Nu (2-59) 
where 
Va 
Y, 
Vie 
[iva] = : (2-60) 
0 
ТЕТІ [0] 
ЕЛЕС m. 3 (2761) 
` [0] | R retor | 
R, 0 0 
N- Mo Re 0 (22,62) 
0 0 R. 
В. ЭК: ЕКЕ 0 55424 0 = RP 
-Ra Rat Ras “Ris 0 0 
Шы > ` : 
~ 0 0 prt Ps Kino Кы 
(2-63) 
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Once again defining a new matrix, the system of equations 


can be expressed as 


d[i] 


IPSI КЕТ (2-64) 
where 
d[L] | | | 
[RR}=[R]+ wee resistance matrix (2-695) 


The equations for the elements of the voltage, 
inductance, and effective resistance matrices are given in 
Appendix A. 

The system of equations describing the electrical 
performance of three-phase induction motor has been 
ed. The elements of the inductance and effective 
resistance matrices are described in terms of the machine 
geometry (airgap active length and width, slot distribu- 
tion, rotor radius, etc.) and electrical parameters 
(number of turns for a stator phase winding, rotor bar 
resistance, and leakage inductances). Thus, in order to 
solve the system of equations for stator phase and rotor 


loop currents, these quantities must be known. 
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2.3 Relation Between System and Single-phase Model 


The object of this section is to derive a relationship 
between the standard single-phase, equivalent-circuit model 
of an induction motor (Figure 2-6) and the system of 
equations developed in the preceding section. From this 
relationship (equations 2-100 to 2-106 of this section), 
the values of the parameters required to solve the system 
of equations can be calculated using the single-phase 


circuit model parameters. 





Figure 2-6. Single-phase equivalent-circuit model for an 
induction machine. 


The parameter values of the single-phase circuit model 


can be obtained from the results of a no-load test, a 


locked-rotor test, and measurements of the DC resistances 
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of the stator phase windings. Section 9-6 of reference [12] 
provides а detailed discussion of the standard test 
procedure and calculations required to obtain the 
equivalent single-phase | circuit parameters for a 
three-phase induction motor. 

In order to derive a relationship between the system of 
equations and a single-phase equivalent circuit model, some 
basic assumptions are necessary. First, the rotor cage is 
symmetric. Thus, each rotor bar is assumed to have the same 
resistance (Ra) and leakage inductance (Шү) Second, 
balanced currents are assumed in both the stator phase 


Windings and rotor loops and are given by 


i, =I,coswt (2-66) 

2n 

ie lucos [i t- 2) (267) 

i t les (o t 21) (2-68) 
= 2 

i=l cos w, teo- RZ) (2-69) 
Nas 


where é is the phase angle of rotor loop 1 current with 
respect to the stator reference frame at t-0. 

Using these assumptions and the system of equations, 
the flux linkages for stator phase a and rotor loop 1 are 


given by the following relations 
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The terminal voltage relations for stator phase a and. rotor 


loop 1 are 


a 
у Зи dt 





+ R¿l¿coswt (22725 


dA, E Oi пр 
0= ——+4R,sin” — |l,cos[w,t+4) (2-73) 
dt RUE 


Defining the variables in complex form according to 


A,= Re(A,exp(jwt)) | (2-74) 
i,= Re(T,exp(jwt)} (2-75) 
v,= Re(V,exp(jwt)) (2-76) 
X, * Re(A,, exp(jo,t]] Que v 
i, * Re{T,exp{jw,t)} (2719) 


and also defining a new rotor loop current in order to 


eliminate phase shifts 
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the terminal voltage relations become 











ESL M Napl 
| Es Lu Jiet jo тмд, (2-80) 
3MI, Nasl пр 
От у +] re о ш 
үш --5 (ҚЫ) ааа (ШЕ), 
АК». пр 
x sie (£5 Yr, (2-81) 


Now, in order to transform the rotor loop equation into 
an equivalent stator phase, a stator-equivalent rotor 
current, lh, is defined. This  stator-equivalent rotor 
current is defined to be the magnitude of the current 
required in the stator windings (under balanced three-phase 
conditions) to produce the same space fundamental component 
of airgap flux as is produced by rotor loop currents of 
magnitude ],,. This corresponds to the equivalent rotor 
current I, of the single-phase, equivalent-circuit of Figure 
2-6. 

The flux linked by stator phase a due to the space 
fundamental component of airgap flux produced by balanced 


rotor loop currents of magnitude lar is 


М Мртав 


Ag 2 


(2-82) 
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The flux linked by stator phase a due to the space 
fundamental component of airgap flux produced by balanced 


stator phase currents of magnitude lI, is 


ST. 
ae (2-83) 





Since the airgap flux must be the same for these currents, 
the following relation between a rotor loop current and the 


equivalent stator phase current is obtained 





3Ls 
= - 84 
Tar (5 =r, (2 ) 


Using this relation to eliminate 1, equations 2-80 and 


2-81 become 








ARI (2-82) 

















[ om | (2-86) 


Multiplying each term in equation 2-81 by the factor 


L;7/M the following result is obtained 
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The voltage relations for the single-phase, equiv- 
alent-circuit model shown in Figure 2-6 are given by 


equations 2-88 and 2-89 below 


Pewee, о о Ба О (2-88) 


R 
О ја Ly T, * 0 (L5 * LT, * —T, (2-89) 


By comparing the terms in equations 2-85 and 2-87 to 
the terms in equations 2-88 and 2-89, the following 


relations are obtained 
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4R m or 
R.= 2sin?( Z2) = 
23 S Nap MON (2-94) 


From the definitions of L,, la, and M given in the 
previous sections, the relationship between the sin- 
gle-phase, equivalent-circuit model parameters and the 


actual machine geometry and electrical parameters is 
R =R, (2-95) 
L, = L... (2-96) 


3 A ,R N?R,d 
pe На s g 


2 
L и < ШЕ -1 (e (2-98) 
NEC sp. ЖЕЛЕЗ ANa J” 
R ee. R 2-99 
2 АШЫ rb ( ) 


From the equations for Liz, і, апа R2, one cannot 
explicitly determine the resistance and leakage inductance 
of a rotor bar. To determine these values the number of 


rotor bars and either the geometry parameters (R,, d, and 
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g) or the number of stator winding turns are required. The 
number of rotor bars can be determined from visual 
inspection or possibly from test methods discussed by 
Hargis et al. [8]. For example, one could use stator frame 
vibration frequencies or rotor slot product frequencies 
present in the stator phase currents. However, determining 
these other parameters will require disassembling the 
machine, removing the rotor, and measuring the rotor 
diameter and length, stator diameter, and number of turns. 
Obviously, this is undesirable. 

Because there is no requirement to determine rotor 
currents explicitly (de. they can't be measured 
directly), it is quite sufficient to deal with equivalent 
rotor currents as seen from the stator windings. This will 
simply result in rotor currents which are scaled by a 
stator-rotor turns ratio and has the effect of rendering 
the choice of N, totally arbitrary. By arbitrarily choosing 
a value for the number stator phase winding turns, the 
rotor loop currents will be scaled by the ratio of the 
actual number of turns to the arbitrary number of turns 
while not affecting in any way the value of stator currents 
or voltages predicted by the model. Although this may not 
be entirely obvious, a simple example will illustrate the 


point. 
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ee 


EXAMPLE: 
Given: A 3-phase, Y-connected, 100 V (line-to-neutral), 60 
Hz, 2-pole induction motor with а 10-bar squirrel-cage 


rotor has the following single-phase, equivalent-circuit 


parameters 


R= 1.00 L,=0.005H L¡,= 0.05H 
L,=0.005H R,-0.50 


The machine is operating at a rotor slip of 0.1 


Equivalent-circuit model solution: 
From equations 2-88 and 2-89, the magnitude of the stator 


phase and equivalent rotor currents calculated are 


1,214.96] A 1,= 13.222 А 


System of equations solution: 


a. Assuming N,-1, the parameters of the system of equations 


are (calculated using equations 2-90 to 2-94) 


В,=1.00 L,,=0.005H L¿= 0.0333H M=0.0131H 
R,,= 0.6760 Lin = 0.0045H L, = 0.024H 
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The magnitude of the stator phase and rotor loop currents 


calculated from equations 2-80 and 2-81 are 


l.=14.961 A lan" 10.082 A 


Using equation 2-84, the stator-equivalent rotor current is 


14-13.222 А 


b. Assuming N,= 10, the parameters of the system of 


equations are (calculated using equations 2-85 to 2-89) 


R,-1.00 Li," 0.005H L,-0.0333H M -0.00131H 
R,,= 0.006760 Li, 7 0.000045H L,=0.00024H 


The magnitude of the stator phase and rotor loop currents 


calculated from equations 2-80 and 2-81 are 


[= 14.961 А [ла = 100.82 A 


Using equation 2-84, the stator-equivalent rotor current is 


[„= 13.222 А 
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The results of (a) and (b) above show that the value of 
the stator current calculated using the system of equations 
is independent of the number of stator turns arbitrarily 
chosen. Although the rotor loop current calculated is 
proportional to the number of stator turns, the 
stator-equivalent rotor current is also independent of the 
number of stator turns arbitrarily chosen. Thus, the number 
of stator turns can be chosen arbitrarily without affecting 
the values of the stator currents calculated using the 
system of equations developed. 

Finally, assuming the number of rotor bars is known and 
arbitrarily setting the number of stator turns to one, the 
parameters required for the system of equations can be 


calculated from the single-phase, equivalent-circuit model 


parameters. 
ee (2-100) 
DeL, (2-101) 
„(8 ја (2- 102) 
ЕЗБЕСІН (2-103) 
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As a final check of the relationship between the system 
of equations and the single-phase, equivalent-circuit 
model, a comparison of the power transferred across the 


airgap will be made. The results must be identical. 


The airgap power for the system of equations is given 


by | ü 


К 


Р атлы р] и (2- 107) 


ier n = bin 1-1) 


Get (2-108) 


- 21, sin P \sin(w,t+0- Nu 


RB 
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The airgap power for the single-phase circuit model is 





TOES (2-398 
(ep) | (2-112) 
ОЗЕ (2-114) 


The airgap power is identical for both the circuit model 
and system of equations. Thus, the relationship between the 


system of equations and equivalent-circuit model is valid. 
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2.4 Summary 


The system of equations describing the electrical 
performance of a three-phase induction motor has been 
developed using stator phase and rotor loop currents as 
state variables. Each stator phase and rotor loop is 
described in terms of resistances and inductances. The 
stator phase voltages are the driving functions for the 
system. The assumptions used to develop the equations 
impose limitations on the ability of the system to 
accurately model any given induction motor. However, each 
assumption can be relaxed and the corresponding system of 
equations can be developed using the method presented. 

A relationship between the standard en 
equivalent-circuit model for an induction motor and the 
system of equations has been derived (equations 2-100 to 
2-106). With this relationship, the electrical and 
geometrical parameters needed to solve the system of 
equations can be easily calculated from the single-phase 
circuit model values. This relationship is useful in that 
the solution to the system of equations can be obtained 
without disassembling the motor, removing the rotor, and 
physically measuring these parameters. Thus, by performing 
the standard no-load, locked-rotor, and DC tests and 
knowing the number of rotor bars, the system of equations 
can be used to simulate any three-phase induction motor, 


subject to the limitations resulting from the assumptions. 
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CHAPTER 3 


COMPUTER SIMULATION OF SYSTEM 


3.1 Requirement for a Simulation Routine 


Due to the number and complexity of the system of 
first-order differential equations developed in Chapter 2, 
a numerical integration routine must be used to solve for 
the time-varying stator phase and rotor loop currents. In 
general, numerical integration methods are used to solve 
systems of equations which are expressed in the following 


form 


d[x(1)] 


de CA E(t) d+ (BI (3-1) 


There are a variety of numerical integration procedures 
and computer simulation programs available which can be 
used to solve the above system of equations. The most 
widely used integration procedures are Euler’s, Milne’s, 
Runge-Kutta, and Hamming’s methods. These methods, which 
are described in references [13,14,15,16,17], are typically 
used to solve initial-value problems. Initial-value 
problems are problems in which the values of the dependent 
variables and necessary derivatives are known at the point 
the integration begins. These various procedures can be 
easily implemented on computers. Several "canned" computer 


programs such as SIMNON, a nonlinear system simulation 
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program developed by the Department of Automatic Control at 
Lund Institute of Technology in Lund, Sweden, and 
PRO-MATLAB, a matrix computation software package by The 
MathWorks, Inc., Sherborn, MA, either use these procedures 
inherently or can be programmed readily to use these 
procedures to solve the system of equations given above. 
The system of equations describing the electrical 


performance of an induction motor is expressed in the form 


d[i(t)]. 


[LI 7 -ERRCOJEQO]* [v 0] (3-2) 


In order to use a numerical integration procedure to solve 
for the time-varying stator phase and rotor loop currents, 
the system of equations must be expressed in the form shown 
in equation 3-1. Comparison of equations 3-2 and 3-3 shows 
that ie inductance matrix in equation 3-2 must be 
inverted. Since some of the elements in this matrix are 
time-dependent, the matrix must be inverted for each value 
of time. This inversion is always possible unless the 
leakage flux is defined to be zero, in which case the 
inverse of the inductance matrix does not exist. With the 
inductance matrix inverted, the system of equations can be 


expressed in the required form i.e., 
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As a result of the inductance matrix inversion 
requirement, SIMNON could not be used to solve the system 
of equations. On the other hand, PRO-MATLAB could be 
programmed to invert the inductance matrix and numerically 
integrate the system of equations. For an induction motor 
with a 45-bar rotor, using a fourth-order Runge-Kutta 
integration procedure, a total CPU run time of five hours 
is required on the LEES Vax-750 to simulate 100 time steps 
using PRO-MATLAB. The major cause of this long run time is 
the inductance matrix inversion. 

In order to reduce the CPU run time required to 
simulate the system, a FORTRAN program was developed. Using 
this program on the LEES Vax-750 for the 45-bar rotor, 100 
time sels. can be simulated in five minutes (approximately 
60 times faster than PRO-MATLAB). In addition to the 
FORTRAN program, which solves the system of equations for 
the time-varying stator phase and rotor loop currents, a 
PRO-MATLAB fast Fourier transform routine is used to 


determine the frequency components of these currents. 


3.1.1 FORTRAN Simulation Program 


From the preceding discussion, the computer simulation 
program must be capable of performing two major steps. 


First, for each value of time, the inductance matrix must 
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be inverted. Second, the resulting system of equations, 
namely equation 3-3 above, must be numerically integrated 
to determine the time-varying values of the stator phase 
and rotor loop currents. Cholesky’s method (also known as 
Crout’s method) with partial pivoting is used to invert 
the inductance matrix. This method was selected over the 
more common Gauss elimination and Gauss-Jordan elimination 
methods since it is more economical with regard to 
computer time and memory allocation [13]. Partial pivoting 
is used to improve the accuracy of the matrix inversion. 
Section 3-4 of reference [13] describes the theory of 
Cholesky’s method. A fixed-step, fourth-order Runge-Kutta 


procedure is used to numerically integrate the system of 


equations. This procedure was selected for several 
reasons. The Runge-Kutta method is a single-step, 
self-starting procedure (i.e., for each integration time 


step, is: the previous value of the dependent variables 
is required). This procedure is relatively easy to program 
and has a small per-step error (on the order of the step 
size to the fifth power). One disadvantage of this method 
is the requirement to perform four function evaluations 
per step. Section 6-5 of reference [13] discusses the 
theory, implementation, and error analysis of the 
fourth-order Runge-Kutta integration procedure. 

The FORTRAN simulation program is divided into six 


parts, a main program and five subroutines. Subject to the 
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assumptions discussed in Chapter 2, the program simulates 
the electrical operation of an induction motor (i.e., 
solves the system of equations for the time-varying stator 
phase and rotor loop currents) with either no broken rotor 
bars or one broken rotor bar. A flowchart for the program 
is shown in Figure 3-1. The FORTRAN source code and sample 
input and output data files are included in section B-1 of 
Appendix B. The following paragraphs describe the purpose 


of each part of the simulation progran. 


Main program: The main program performs several 
functions. The primary function is to act as a buffer 
and pass data between the subroutines. In addition, 
the main program keeps track of the simulation time, 
calculates the value of the current vector at the end 
of each time step, and writes the results to an 


external file. 


Subroutine INPUT: This subroutine is called once by 
the main program. It performs two functions. First, it 
reads the data contained in an external input file. 
The external input file provides a description of the 
machine being simulated. It includes the single-phase, 
equivalent-circuit model values, number of rotor bars, 
and number of pole-pairs for the machine, the rotor 


slip, and the supply voltage and frequency. The 
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Gun 


declare variables 


input parameters 
INPUT 


initialize variables 
2003720, tjme-0 









calculate system matrices 
for time 
MOTOR 


start simulation loop 


calculate [v]-[RR}[i] 
DRIVE 


first R-K factor 


RK1=d/dt[i] 
CHLSKY 


time=t+inc/2 


calculate intermediate current 
[EST]z[i]*REKl*inc/2 
CURRENT 








calculate system matrices 
for time 
MOTOR 






calculate [v]-[RR][EST] 
DRIVE 


second R-K factor 
RK2=d/dt[i] 


CHLSKY 





2) 


Figure 3-1. FORTRAN simulation program flowchart. 
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calculate intermediate current 
[EST]=[i]+RK2*inc/2 
CURRENT 


calculate [v)]-[RR][EST)] 
DRIVE 


third R-K factor 


RK3zd/dt[i] 
CHLSKY 


calculate intermediate current 
[ESTJ=[i)J+RK3tinc 
CURRENT 














calculate system matrices 
for time 
MOTOR 


calculate [v]-[RRJ[EST] 
DRIVE 


fourth R-K factor 
RK4zd/dt[i] 
CHLSEY 


calculate new current values 
and increment simulation time 

[iJ=[iJ+(inc/6) (RK1+2*RK2+2*RK3+RK4 ) 
t=t+inc š 





write t,[i] to output file 


<--- true 





ae 


false 


Figure 3-1.(cont.) FORTRAN simulation program flowchart. 
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simulation time parameters (step-size and stop time) 
and a broken bar flag are also included in the input 
file. Second, using the relations derived in section 
2.3, the coefficients for the inductance and effective 
resistance matrix elements are calculated and passed 


to the main progran. 


Subroutine MOTOR: This subroutine is called for each 
value of time. It performs two functions; Date 
calculates the elements of the voltage, inductance, 
and effective resistance matrices, and it reduces the 
system of equations. For the voltage vector, a 
balanced three-phase supply is assumed. The inductance 
and effective resistance elements are calculated using 


the relations given in Appendix A. For both a "no 


simulation and a one broken bar 


broken bar 
ación, the system of equations can be reduced 
from Nr+3 to №+2 since the sum of the rotor loop 
currents equals zero. This condition is imposed on the 
system of equations by eliminating rotor loop Nis 
current. Column N,+3 is subtracted from columns 4 
through N,,+2 and column N„+3 is deleted from the 
inductance and effective resistance matrices. Next, 
row Na+3 is deleted from each matrix to complete the 


" 


reduction. In addition, for а one broken bar 


" 


simulation (assume rotor bar N»s-1 is broken), the 
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system of equations can be reduced again by setting 
rotor loop Nas-1 current equal to rotor loop Nas-2 
current. For this condition, column N,,*2 is added to 
column Natl and column N,+2 is deleted from the 
inductance and effective resistance matrices. Fara- 
day's law is satisfied by adding row Nast2 to row Natl 
and deleting row N,,*2 from each matrix to complete the 


reduction. 


Subroutine CURRENT: This subroutine is used to 
calculate the intermediate values of the current 


vector for the Runge-Kutta method. 


at) 


HSA: 


At (3-4) 


where 


[EST]= intermediate value of current vector 


1] = initial current vector 


Subroutine DRIVE: The purpose of this subroutine is to 
calculate the "driving function" for each value of 


time and current, i.e., 


driving function- [V]- [RR ][1] (325) 
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Subroutine  CHLSKY: This subroutine performs the 
inductance matrix inversion to calculate the deriva- 
tive of the current vector for each time value. For 
the fourth-order Runge-Kutta procedure, this calcula- 


tion is performed four times per time step. 


The program can be modified, if desired, to simulate a 
variety of operating/fault conditions. For example, an 
unbalanced voltage supply can be simulated by simply 
changing the appropriate statements in subroutine MOTOR. 
Similarly, the effects of an increase or decrease in the 
impedance of one or more rotor bars can be simulated by 
modifying the appropriate statements for the elements of 
the inductance and effective resistance matrix in 
subroutine MOTOR. In addition, multiple broken rotor bars 
can be simulated by repeating the system reduction scheme 


in subroutine MOTOR for each broken bar. 


- 


3.1.2 Determining the Time Step 


An appropriate time step for any simulation can be 
determined using a "trial and error" procedure. Basically, 
a series of simulations are run with increasing values of 
time steps and the output results are compared. The 
largest time step that gives "good" results (i.e., same 
output values as those using a smaller time step) is an 


appropriate time step. 
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Another method which provides an estimate for an 
appropriate time step and also provides an indication of 
the system stability and of the time required for the 
solution to reach a steady-state value is to determine the 
eigenvalues of the system when the rotor speed is equal to 
zero. From equation 3-3, the eigenvalues of interest are 
those of the matrix -[L]-![RR]. For a rotor speed of zero 
(s=1.0), this matrix is time-invariant and the system is 
stable if all the eigenvalues are negative. The inverse of 
the magnitude of the largest eigenvalue represents the 
“fastest” transient time for the system and thus provides 
an estimate for the largest time step. Similarly, the 
inverse of magnitude of the smallest eigenvalue represents 
the "slowest" transient time for the system and thus 
provides an estimate of the time required for the solution 
to reach steady-state. PRO-MATLAB, which has a built-in 
SS value function, can be used to determine these 
eigenvalues. A PRO-MATLAB routine which calculates the 
system eigenvalues for s=1.0 and the eigenvalues for a 
3-bar rotor and а 45-bar rotor induction motor are 


included in section B-2 of Appendix B. 


3.1.3 Running the Simulation Program 


In order to start the simulation program, the stator 
phase and rotor loop currents must be initialized. Since 
the steady-state values of the stator phase and rotor loop 


currents for a "one broken bar" case are not known, the 
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currents are set to zero initially and the simulation 
program runs until steady-state conditions exist. For the 
3-bar rotor and 45-bar rotor induction motors simulated, 
steady-state conditions were reached after approximately 


0.2 seconds of simulation time. 


3.1.4 Processing the Simulation Output 


In order to determine the harmonic content of the 
time-varying stator phase and rotor loop currents, a 
transformation from the time to the frequency domain is 
necessary. This transformation is accomplished by taking 
the fast Fourier transform (FFT) of the time-varying 
stator phase and rotor loop currents. For a broken rotor 
bar, the stator phase currents will contain a harmonic 
component at a frequency of (1-2s)f (refer to section 
1.3). Most induction machines are designed to operate at 
small slip values on the order of 0.02-0.05. Thus, the 
frequency of this harmonic component is very near the 
fundamental frequency (0.9f to 0.96f). To distinguish this 
component from the fundamental frequency component using 
an FFT routine, a frequency resolution on the order of 
0.02f is required. From reference [18], the bin width ds 


frequency resolution of a FFT for N samples is given by 


Af - = 222) 
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where 


Af - frequency resolution 
f. "sample frequency 


N=number of data samples 


In addition to selecting the proper frequency 
resolution, a weighting function, or window, can be 
applied to the sampled data to aid in distinguishing 
frequency components which are relatively close. Refer- 
ences [18,19] provide a detailed description of the fast 
Fourier transform and common window functions used, their 
properties, and their advantages and disadvantages. 

PRO-MATLAB, which has a built-in FFT function, is used 
to transform the stator phase and rotor loop currents from 
the time domain to the frequency domain. A frequency 
resolution of 0.3 Hz and a second-order Hanning window 
[18] are used to enhance the detection of the (1-2s)f 
component of the stator phase current. A sample frequency 
of 153.6 Hz and a sample size of 512 data points satisfy 
the 0.3 Hz frequency resolution requirement and PRO-MAT- 
LAB’s requirement for a sample size which is an integer 
power of two. The PRO-MATLAB FFT routine and a sample 


output file is included in section B-3 of Appendix B. 
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The relation between the magnitude of the FFT using a 
second-order Hanning window and the magnitude of the 
corresponding periodic signal in the time domain is given 


by [18] 


N 
IFCoo)| * 7I £ CO] (3-7) 


where 


\F(w)|= magnitude of FFT at w 


f(t)-* |f(t)] cos(ot* 0) 


For the PRO-MATLAB FFT routine, the magnitude of the 
FFT is converted to a dB reference scale (for graphing 


purposes) using the relation 


dB= 20log|F(w)| (3-8) 


3.2 "Hand" Verification of Simulation 


For initial validation and "debugging" of the FORTRAN 
simulation program and FFT routine, a two-pole, 3-bar rotor 
induction motor was chosen. An exact solution for the 


stator phase and rotor loop currents can be calculated by 
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hand for this simple machine with either no broken rotor 
bars or one broken rotor bar. The exact solutions can be 


calculated using the method of undetermined coefficients. 
No broken rotor bars: The system of equations for this 


case can be reduced from six to four by imposing the 


following constraints 


jeep (3-9) 


iy +#1,2+1 7 0 (3-10) 


Assuming the following form for the stator phase and rotor 


loop currents and substituting into the system of 


equations, the unknown coefficients (A,B,C,D) can be 
determined. 

1,7 Acoswt+Bsinwt (3-17) 

2n 
1, = a cos 01-22) 8510 wt- E) (3212) 

3 3 
i,, = Ccos(swt)+ Dsin(swt) (3- 13) 

2m 

iam Ccos(swt- 75)» Dsin( sat-22) (3-14) 
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The amplitude of the stator phase and rotor loop currents 


are given by the following relations 


(60 Hz)* li,| * li,| - J A?« B? (3-15) 


Ire [i| - Jil - 4 C? D? (3295) 


ne broken rotor bar: Assuming rotor bar 2 is the 
broken bar, the system of equations can be reduced from six 
to three by imposing the constraints given by equations 3-9 


and 3-10 plus the constraint 
UY CT E (ЕУ) 


Again, assuming the following form for the stator phase and 
rotor loop currents, the unknown “coefficients can be 


determined. 


i," Acoswt+Bsinwt+ 


Ccos((1-2s)w t) * Dsin((1-2s)o t) (3-18) 


2 
і, = Acos( wt- 2") Bsin( w em). 


ces (11 2) 1-22) osis(ai- 2991-57] (3-19) 
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i,, = Ecos(swt)+Fsin(swt) (3 — 20) 


The amplitude of the stator phase currents (both the 60 Hz 
and (1-2s)f components) and the rotor loop currents are 


given by the following relations 


1(60 Hz) = |i,(60Hz)|=|i,(60Hz)]= /A?+B? (3-21) 
I(1-2s)f - [i,(1-2s)f| - [i,(1- 2s)f| - 4 C? D? (3-22) 


геј |- E? F? (3- 23) 


3.2.1 Three-bar Rotor Simulation 


To simulate a 3-bar rotor induction motor, the 
single-phase, equivalent-circuit model parameters from the 
3-HP experimental motor are used as inputs for the FORTRAN 
Simulation program. The parameters used for the 3-bar 
rotor induction motor are listed in Table 3-1. 

Using the method described in the preceding рага- 
graphs, the exact solution of the stator phase and rotor 
loop currents were calculated for various rotor slip 
values. Table 3-2 shows the results from these 
calculations. The results generated from the simulation 


program are shown in Table 3-3. 
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РАКАМЕТЕК 


0.8590 ohms 
0.0046 H 
0.0704 H 
0.0046 H 

0.5612 ohms 

3 
1 


Input Frequency 60 Hz 





Input Phase Voltage 169.71 V 


Table 3-1. 3-bar rotor parameters. | 


The correlation between the exact and simulation 
results is extremely good. For all cases, the error 
between these results is less than 1%. The absolute error 
between the simulation and exact results is shown in Table 
3-4. 

The close correlation between the exact and simulation 
results is demonstrated graphically in Figures 3-2, 3-3, 
3-4, 3-5, 3-6. Based on the results presented above, the 
FORTRAN simulation program and FFT routine can be seen to 
accurately solve the system of equations developed in 


Chapter 2 for the stator phase and rotor loop currents. 
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| NO BROKEN BARS | | NO BROKEN BARS | BARS ONE BROKEN BAR 
SLIP |1(60 На) I(60 Hz) | I(1-2s)f 
Amps Amps Amps Amps Amps 


23.98 20.02 15.42 11553 о 1 


46.71 39.76 41.58 Eu: 


Table 3-2. 3-bar rotor exact results. 


NO BROKEN BARS ONE BROKEN BAR 
SLIP | 1(60 Hz) | ; (60 Hz) | I(1-28)f 
Amps ш Amps Amps Amps 
23.85 19.87 19237 11.49 11.44 


46.71 39.75 41.58 Fr 19.87 


Table 3-3. 3-bar rotor simulation results. 












19.88 


N Rar de exp EE =< 
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E 







! NO BROKEN BARS | ONE BROKEN BAR 
: Ir 


Ir | 





Table 3-4. Error between exact and simulation results. 


For this simple machine, a broken rotor bar 
represents a significant change in the rotor magnetic 
circuit. As a result of this change, the stator phase 
current is also significantly altered. The drastic change 
of the stator phase current (in both the time and 
frequency domain) for a rotor slip of 0.05 is shown in 


figures 3-7, 3-8, 3-9, 3-10. 
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(6082) in AMPS 


A ESR ROTOR STATOR PHASE CURRENT 
SIMULATION vs. EXACT RESULTS 
NO BROKEN BARS 





0.00 0.20 0.40 0.60 0.80 1.00 
ROMORES El 


— simulation values 
-X€- exact values 


Figure 3-2. I(60 Hz) with no broken bars. 
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I(rolor loop) in AMPS 


J BAR ROTOR ROTOR LOOP CURRENT 
SIMULATION vs. EXACT RESULTS 
NO BROKEN BARS 





0.00 0520 0.40 0.60 0.80 1.00 
ROTOR SLIP 


simulation values 
exact values 


Figure 3-3. Ir with no broken bars. 
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(6014) in AMPS 


M ROTOR STATOR PHASE CURRENT 
SIMULATION vs. EXACT RESULTS 
ONE BROKEN BAR 





0.00 0.20 0.40 0.60 0.80 1.00 


Figure 3-4. 


RO TOR Sie 


simulation values 
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I(60 Hz) with one broken bar. 
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I(1-2s)f COMPONENT in AMPS 





Soe KOTOR STATOR PHASE CURRENT 
SIMULATION vs. EXACT RESULTS 
ONE BROKEN BAR 


+ 


0.00 (ЙО 0.40 0.60 0.80 1.00 
ROTOR SLIP 


0 simulation values 
x exact values 


Figure 3-5. I(1-2s)f with one broken bar. 


98 





(го(ог loop) in AMPS 


coo ROTOR ROTOR LOOP CURRENT 
SIMULATION vs. EXACT RESULTS 
ONE BROKEN BAR 





0.00 (2240 0.10 0.60 0.80 1.00 
КОЛОС ЗА 


sımuflazsımn values 


< exact values 


Figure 3-6. Ir with one broken bar. 
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SIMULATION STATOR PHASE CURRENT 


.05 


3—BAR ROTOR slip 


NO BROKEN BARS 





SdNV WU! ENC) 


.167 


1 


2103.3 


1.067 1.100 1 
TIME in SECS 


1.033 


1.000 


Stator phase current with no broken bars. 


Figure 3-7. 
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CURRENT in AMPS 


Fi 


gure 


SIMULATION STATOR PHASE CURRENT 
3-BAR ROTOR slip=.05 
ONE BROKEN BAR 





1.000 1.067 1.133 1.200 1726: 1.333 
PAE sin SECS 


3-8. Stator phase current with one broken bar. 





STATOR PHASE FREQUENCY SPECTRUM 
3 BAR ROTOR Slip 05 
NO BROKEN BARS 


SPECTRAL POWER DENSITY in dB 





4 56 58 60 62 64 66 68 70 
FREQUENCY in Hz 


S0 De 5 


Figure 3-9. Frequency spectrum with no broken bars. 
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ПРИ МИЛИ ОРАО АНИ И О РАСТАО 
3-BAR ROTOR slip=.05 
ONE BROKEN BAR 


POWER DENSITY in dB 


SPECTRAL 





14 об a8 60" 62 64° 66 68 970 
FREQUENT cel 


ANS 


Figure 3-10. Frequency spectrum with one broken bar. 
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3.3 Simulation Results For Experimental Motor 


Now that the simulation program has been initially 
validated, it can be used to predict the results of a real 
induction motor. These predicted results will then be 
compared to actual experimental results to validate the 
methodology presented in Chapter 2 and the simulation 
program. A series of experiments will be conducted using a 
3-phase, 3-HP, 230 V (line-to-line), 60 Hz, 1730 rpm 
"off-the-shelf" induction motor. Due to testing facility 
limitations (explained in Chapter 4), the experiments can 
only be conducted for small values of rotor slip (0 to 
0.04). Although this is a rather small range, it does cover 
the full load operating range of the experimental motor. In 
addition, only one rotor bar will broken. The following 
sections present the simulation results for the experimen- 
tal motor operating at small slip values with and without a 
broken rotor bar. For the "no broken bar" simulations, the 
simulation results are compared to the results calculated 


using the single-phase, equivalent-circuit model for this 


machine. The parameters required to simulate the 
experimental motor are listed in Table 3-5. The 
single-phase, equivalent-circuit model parameters were 


determined using а least squares parameter estimation 


routine [20]. 
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PARAMETER VALUE 


0.8590 ohms 












0.0046 H 
0.0704 H 
0.0046 H 
0.5612 ohms 
45 


2 







Input Frequency 60 Hz 





Input Phase Voltage 169.71 V 


Table 3-5. Experimental motor parameters. 


3.3.1 Case 1: No Broken Rotor Bars 


The simulation results, equivalent-circuit model 
results, and the error between these results are given in 
Table 3-6. As expected, the error between the simulation 
and equivalent-circuit model results is less than 0.5% for 
each rotor slip value. Figure 3-11 shows the small 


difference between these results. 
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SIMULATION EQ. CKT. MODEL * ERROR 


I(60 Hz) I(60 Hz) 





Table 3-6. Simulation results-no broken bars. 
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SIMILATION RESULTS FOR TEST MOTOR 
STATOR PHASE CURRENT 
NO BROKEN BARS 


10 


(60 Пи) in AMPS 





0.00 0.01 0.02 0.03 0.04 
ROTOR SLIP 


simulation results 
equivalent model values 


— 
— 
` 
SS o 


Figure 3-11. Simulation vs. equivalent circuit results. 
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The simulation results for a one broken bar case are 
given in Table 3-7. Plots of the fundamental and (1-2s)f 
component of the stator phase current as a function of 
rotor slip are shown in Figures 3-12 and 3-13. The results 
show that the (1-2s)f component of the stator phase 
current increases approximately linearly with increasing 
rotor slip. Unlike the 3-bar rotor case above, there is 
only a slight change in the fundamental component with 


respect to the no-broken-bar case. This is demonstrated 


graphically in Figure 3-14. 


SLIP I(60 Hz) I(1-2s)f I(1-2s)f/1(60 
| Amps Amps 





Table 3-7. Simulation results-one broken bar. 
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(60 Hz) in AMPS 


SIMWLATION. RESULTS FOR TEST MOTOR 


ШАРЫ PHASE CURRENT 
ОМЕ ВКОКЕМ ВАК 


0.00 0.01 0.02 0.03 0.04 
ROTOR SLIP 


== simulation results 


Figure 3-12. I(60 Hz) vs. slip. 


109 





I(1—2s)f in AMPS 


SIMULATION RESULTS FOR TEST MOTOR 
ПАОЛА А У I 1. 2s fe CORRENT 
ONE BROKEN Bak 
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= 
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0.00 


0.00 0.01 0.02 0.03 0.04 
ROTOR SLUF 


—- simulation results 


Figure 3-13. I(1-2s)f vs. slip. 
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I(60 Hz) in AMPS 


SIMULATION RESULTS FOR TEST MOTOR 
STATOR PHASE CURRENT 
GOOD vs BAD ROTOR 


10 





0.00 0.01 ШЕШЕ 0.03 0.04 
ROTOR SLIP 


—— no broken bars 
-ж- one broken bar 


Figure 3-14. I(60 Hz) no vs. one broken bar. 
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3.4 Summary 


A FORTRAN simulation routine has been developed to 
solve the system of equations describing the electrical 
performance of an induction motor derived in Chapter 2. The 
routine uses a fixed-step, fourth-order Runge-Kutta 
procedure to numerically integrate the system of equations. 
A matrix inversion routine, Cholesky's method with partial 
pivoting, is required to invert the time-varying inductance 
matrix for each value of time. The program can be used to 
simulate an induction motor with either no broken rotor 
bars or one broken rotor bar. In addition to the FORTRAN 
simulation program, two  PRO-MATLAB routines, are also 
employed. The PRO-MATLAB eigenvalue routine is used to aid 
in determining the fixed time step for the Runge-Kutta 
integration. A PRO-MATLAB FFT routine is used to transform 
the output data of the FORTRAN program from the time domain 
to the - frequency domain. Thus the RES frequency 
components of the stator phase and rotor loop currents can 
be determined. Appendix B includes the source code and 
sample input and output files for these routines. 

The FORTRAN simulation program and PRO-MATLAB  FFT 
routine were initially validated using a 3-bar rotor 
induction motor. The error between the simulation results 


and exact solutions for both a no-broken-bar and a 





one-broken-bar case is less than 1%. The simulation program 
accurately solves the system of equations derived in 
Chapter 2 for the stator phase and rotor loop currents. 

The simulation results for the experimental motor have 
been presented and will be compared to actual experimental 
results for final validation of the simulation program in 


Chapter 4. 





СНАРТЕК 4 
INVESTIGATION OF AN INDUCTION MOTOR 


WITH A BROKEN ROTOR BAR 


4.1 Introduction 


A series of experiments was conducted using a set of 
initially identical “off-the-shelf” J-HP  squirrel-cage 
induction motors with various rotors modified so as to 
investigate their performance with and without a broken 
rotor bar. The magnitudes of the fundamental and (1-2s)f 
harmonic components of the stator phase current were 
recorded. Two identical, "as manufactured" rotors were used 
to collect data corresponding to a "good" motor with no 
broken rotor bars. A third rotor, identical to the other 
two with the exception that a single rotor bar was 
deliberately open-circuited, was used to collect data 
corresponding to a "broken rotor bar" motor. А single 
stator I used with the rotors being swapped in and out 
for the various tests. 

The nameplate data for the test motors is given in 
Table 4-1. As stated in the preceding paragraph, three 
rotors were alternately installed into the same stator for 
the experiments. Table 4-2 provides a summary for these 
rotors. For the broken bar rotor, one end of a single rotor 
bar was initially "disconnected" from the end ring by 
removing a small section of the rotor bar at the rotor 


bar-to-end ring interface with a milling machine. After the 





initial tests were completed, the opposite end of the same 
rotor bar was also "disconnected" from the other end ring 
and a series of tests were again conducted using this rotor 


(ROTOR $2.1). 


Westinghouse Life-Long T AC Motor 


3-phase 60 Hz 3-HP 1730 rpm 









230/460 V 9.4/4.7 A 182T frame 1.0 s.f. 





catalog no.: 05-3H4SBFC-SKB motor style: 773B646G41 








Table 4-1. Test motor nameplate data. 


“еее | 
[eee | o mn mn | 
| пета a [sen a ea aee | 
As manufactured; no broken bars. 


Table 4-2. Rotor summary. 















The experimental data collected for each motor serves 
two primary purposes. First, comparison between the 
experimental results and predicted results (presented in 


Chapter 3) will be used to validate the computer simulation 





program. Of interest is whether the simulation program 
accurately solves for the stator phase currents of an 
induction motor with and without a broken rotor bar. In 
addition, this comparison can provide an indication of the 
areas which must be added and/or modified, if required, in 
the simulation program to accurately describe the 
electrical performance of an induction motor with ог 
without a broken rotor bar. Second, comparison of the 
experimental results for a "good" motor and a "broken rotor 
bar" motor can be used to make a quantitative decision 
regarding the feasibility of detecting broken rotor bars by 
monitoring the magnitude of the (1-2s)f harmonic component 


of the stator phase current. 


4.2 Description of Experimental Facility 


The facilities used to conduct the experiments consists 
of four major components; a 208 V rms (line-to-line) 
3-phase oe supply, a dynamometer, a test motor, and a 
Signal analyzer. A MAGTROL, Inc. model HD800-8 dynamometer 
is used to provide the load necessary to operate the test 
motor at various constant speeds (or rotor slips). Due to 
the kinetic power dissipation limit of the HD800-8, only 
small values of rotor slip (0-0.04) can be safely achieved. 
However, this range of rotor slips does cover the full load 
operating range of the test motors. For a rotor slip of 


0.04, the test motor can be operated for approximately 10 





minutes every two hours without exceeding the rating of the 
dynamometer. The control section for the HD800-8 provides 
digital read-outs for the following parameters: 

stator phase current [A rms] 


stator phase voltage [V rms] 


stator phase power [W] 
shaft speed [rpm] 
shaft power [W] 
shaft torque [N-m] 


A Hewlett-Packard model HP-3561A dynamic Signal 
analyzer is used to fast Fourier transform the stator phase 
current measurements and display the resulting frequency 
spectrum. The HP-3561A is a single-channel, multi-function 
Signal analyzer. The analyzer was set-up to calculate and 
display the stator phase current frequency spectrum over a 
20 Hz span centered at 60 Hz. For this setting, the 
frequency resolution of the analyzer is 0.05 Hz and a 
sample-time record length of 20 seconds is required to 
process the stator-current signal, Like the simulation 
program, a second-order Hanning window function is used 
to aid in distinguishing among frequency components which 
are relatively close. In order to increase the sig- 
nal-to-noise ratio of the input stator-current signal, a 
time averaging feature is used. Using this feature, five 


successive time records are averaged on a point-by-point 





basis. Thus, for each test the analyzer requires 
approximately two minutes of stator phase current signals. 
The output units of the HP-3561A are volts rms squared. For 
the “typical” frequency spectrum plots included in this 
chapter (Figures 4-4, 4-8, 4-15, and 4-19), the following 
relation can be used to convert the analyzer output to amps 


rms 


A rms=2.0y(V rms)? (4-1) 


where 


(V rms)/* analyzer output reading 


A major concern throughout the experimental portion of 
this research is the issue of data reproducibility. In an 
effort to increase the confidence level of the experimental 
data seis a, several tests were conducted PE each rotor 
at each rotor slip value. For each test the magnitudes of 
the (1-2s)f and fundamental components of the stator phase 
currents were recorded. These values were then averaged and 
the standard deviation was determined. The tests were 
continually repeated until the measured values were 
consistently within one standard deviation of the mean 
value. In addition to this, the averaging technique, number 


of averages, frequency span and resolution, etc. of the 
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HP-3561A were varied to ensure that the choice of analyzer 
settings was not introducing any error into the 


measurements. 


— — SSS, U 


The following sections provide a summary of the 
experimental data collected and corresponding analysis for 
each rotor. Appendix C contains tables of the measured data 
(the fundamental and (1-2s)f components of the stator phase 


current) for each rotor tested. 


4.3.1 Experimental Results for a Motor Without a Broken 
Rotor Bar(ROTOR #1) 


The averaged voltage and current data for ROTOR #1 is 
shown in Table 4-3. This is the baseline data. In addition 
to the fundamental and (1-2s)f components of the stator 
phase current, three other components were apparent, at 
frequencies (1+2s)f, (1-4s)f, and (1+4s)f. In general, the 
(1-4s)f and (1+4s)f components, when detected, were 
usually on the order of one-half the magnitude of the 
(1-2s)f component. As shown in Appendix D, these 
components could be produced by the third and fifth 
time-harmonics of the stator phase current. The (1+2s)f 
component is typically on the same order as the (1-2s)f 
component and is generally attributed to torque pulsations 


and the resulting speed oscillations at twice-the-slip 
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frequency [7,8,9,10,11]. In addition, this component could 
also be produced by the third time-harmonic of the stator 
phase current 85 shown in Appendix D. Table 4-4 lists the 
magnitudes of these components for в "typical" test at 
each each rotor slip value. 

For all values of rotor slip, the non-fundamental 
harmonic components detected were less than 0.05% of the 
fundamental component of the stator phase current. From 
Table 4-3, the variation in the magnitude of the (1-2s)f 
component is shown to be on the order of 25-30% of the 
mean value for this component. Figures 4-1 and 4-2 show 
the averaged data and raw data for the (1-2s)f component. 
Figure 4-3 shows the averaged data for the fundamental 


component. Figure 4-4 shows the "typical" frequency 


Spectrum for various rotor slip values. 







I(1-2s)f 


SLIP VOLTAGE STD. DEV. STD. DEV. 
A rms A rms 


V rms 
пи | зна | км | ом pem em 
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ale Le L o om [= 


Table 4-3. ROTOR #1 averaged data. 
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Table 4-4. ROTOR £1 "typical" harmonic frequency data. 


121 
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Figure 4-1. ROTOR #1 I(1-2s)f averaged results. 
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Figure 4-2. ROTOR #1 I(1-2s)f "raw" data. 
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ROTOR #1 I(60 Hz) averaged results. 
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Figure 4-4. ROTOR £1 "typical" frequency spectrum. 
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Bar (one end open--ROTOR #2) 


The averaged data for ROTOR #2 is shown in Table 4-5. 
Like the results obtained for ROTOR #1, the components at 
frequencies (1+2s)f, (1-4s)f, and (1+4s)f were detected in 
addition to the fundamental and (1-2s)f stator phase 
current harmonics. The (1-4s)f and (1+48)f components, 
when detected, were usually on the order of one-half the 
magnitude of the (1-2s)f component. The (1+2s)f component 
was typically on the same order as the (1-2s)f component. 
Table 4-6 lists the magnitudes of these components for a 
“typical” test at each rotor slip value. 

For all values of rotor slip, the non-fundamental 
harmonic components detected were less than 0.05% of the 
fundamental component of the stator phase current. From 
Table 4-5, the variation in the magnitude of the (1-2s)f 
component is shown to be on the order of 25-30% of the 
mean value for this component. Figures 4-5 and 4-6 show 
the averaged data and raw data for the (1-2s)f component. 
Figure 4-7 shows the averaged data for the fundamental 
component. Figure 4-8 shows the "typical" frequency 


spectrum for various rotor slip values. 


126 





I(60 Hz) I(1-2s)f 
SLIP VOLTAGE А š 
V rms 
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Table 4-5. ROTOR #2 averaged data. 
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Table 4-6. ROTOR $2 "typical" harmonic frequency data. 
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EXPERIMENTAL RESULTS 
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Figure 4-5. ROTOR #2 1(1-2s)f averaged results. 
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Figure 4-6. ROTOR #2 I(1-2s)f "raw" data. 
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Figure 4-7. ROTOR $2 I(60 Hz) averaged results. 
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Figure 4-8. ROTOR #2 "typical" frequency spectrum, 
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4.3.3 Discussion of Simulation and Experimental Results 


Tables 4-7 and 4-8 show the simulation results from 
Chapter 3 and experimental results of the fundamental and 
(1-2s)f harmonic components of the stator phase currents 
for a motor with and without a broken rotor bar. These 
tables show that the simulation program accurately solves 
for the fundamental component of the stator phase current 
for a motor with and without a broken rotor bar. This was 
expected since the input parameters for the simulation 
program are determined using the single-phase, equiv- 
alent-circuit model values for the test motors. However, 
for the (1-2s)f component, the simulation results and 
experimental results are significantly different. These 
unexpected differences raise two questions which must be 
answered. 

First, why is the value of the (1-2s)f component of 
the stator phase current non-zero in a motor without a 
broken rotor bar? This question can be answered by 


considering the effect of manufacturing asymmetries on the 


stator phase current. From the analysis conducted by 
Kliman et al. [7], it was shown that manufacturing 
asymmetries (i.e., rotor out of round, rotating eccentric- 


ity, or non-uniform magnetic orientation of the rotor 
laminations) can also create airgap flux anomalies with a 


fundamental component on the same order as that of a 
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broken bar. This flux anomaly induces currents at a 
frequency of (1-2s)f in the stator phase windings. Thus, 
manufacturing asymmetries, which are present to some 
degree in all induction motors, can result in a non-zero 


(1-2s)f harmonic component of the stator phase current. 


1(60 Hz) AMPS rms 


"GOOD" MOTOR "BROKEN BAR" MOTOR 
ROTOR #1 SIMULATION ROTOR #2 SIMULATION | 


4.76 4.79 


.95 





Table 4-7. Simulation and experimental results 1(60 Hz). 


I(1-2s)f AMPS rms 


o = = 
"GOOD" MOTOR “BROKEN BAR" MOTOR 
| ROTOR #1 SIMULATION ROTOR #2 SIMULATION | 
3 


0.001 


0.0018 
0.0023 





Table 4-8. Simulation and experimental results I(1-2s)f. 
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The second question which must be answered is why are 
the simulation results for the (1-2s)f component on the 
order of 40 to 75 times greater than those measured for a 
motor with a broken rotor bar. One possible mechanism is 
that the bar was not completely broken. This could be the 
case in the present experiments if the rotor bar was not 
completely disconnected from the end ring during milling 
so that path for current to flow still exists. In order to 
demonstrate how this mechanism would result in a reduced 
magnitude of the (1-2s)f component, the simulation program 
was used to calculate the stator phase currents for a 
motor with a rotor bar that has a higher series impedance 
than the other rotor bars. For a rotor slip of 0.04, the 
impedance (both the rotor bar resistance and leakage 
inductance) of one rotor bar was varied (relative to the 
impedance of a normal rotor bar) and the fundamental and 
(1-2s)f components of the stator phase current was 
calculated. These results are shown in Table 4-9 and 
graphically in Figures 4-9 and 4-10. As expected, the 
results show that the fundamental component of the stator 
phase current remains approximately constant while the 
Magnitude of the (1-2s)f component increases as the rotor 
bar impedance is increased. In addition, these results 
indicate that only a 2% increase in the bar impedance is 
required for the simulation program to predict the (1-2s)f 


values measured for ROTOR #2. Certainly, the bar impedance 
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would increase by a substantially larger factor than 2% if 
a section of the rotor bar was removed (i.e., the broken 
bar on ROTOR #2). 

After removing ROTOR #2 and visually inspecting the 
broken rotor bar, it was concluded that the rotor bar was 
completely cut away from the end ring. Thus it appears 
certain that current is not flowing into the broken bar 


across the break from the end ring. 





RATIO OF BAR IMPEDANCE 1(60 Hz) I(1-28)f | 
TO ALL OTHER BARS AMPS rms AMPS rms | 
i | Р | ( ma и и 
















infinity 0.1758 
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NOTE: all values for rotor slip=0.04 


Table 4-9. Bar impedance variation results. 
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Figure 4-9. Bar impedance variation I(1-2s)f. 
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Figure 4-10. Bar impedance variation I(60 Hz). 
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Another mechanism that could result in a current 
flowing through a broken rotor bar is that of inter-bar 
currents. Inter-bar currents are currents which flow 
between a rotor bar and adjacent rotor bars through the 
rotor laminations. The presence of such currents could be 
expected to "short circuit" the effect of a break in a 
rotor bar, hence resulting in a current flow in the broken 
bar that would not be predicted by the simulation program. 

The existence of large inter-bar currents in 
three-phase squirrel-cage induction motors with broken 
rotor bars has been analyzed and experimentally verified 
by Kerszenbaum and Landy [9]. They have shown that the 
Magnitude of the inter-bar current is a function of the 
rotor-bar impedance and the distributed inter-bar 
impedance. The inter-bar impedance is the sum of the 
contact impedance of the rotor bar-to-rotor lamination 
interface and the impedance of the rotor core (i.e., 
laminations) between two adjacent rotor bars. As the ratio 
of the rotor-bar impedance to the КЕ impedance 
increases, the magnitude of the current flowing into a 
broken rotor bar approaches the magnitude of the current 
which would flow into a "healthy" (i.e., not broken) rotor 
bar for the same operating conditions. This is shown in 
Figure 4-9 where 
= current flowing in a broken bar 


m 


m = current flowing in a "healthy" bar 
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Figure 4-11. Ratio of current in a broken bar to a 
"healthy" bar [9]. 


Bar (both ends open--ROTOR #2.1) 


In order to investigate the hypothesis that inter-bar 
currents were masking the effects of the break in the rotor 
bar during experiments with ROTOR #2, the opposite end of 
the broken rotor bar on ROTOR #2 was also cut away from the 
end ring, resulting in the rotor known as ROTOR #2.1. This 
was expected to result in a further decrease in the 
Magnitude of the current flowing through the broken rotor 
bar thus, increasing the resulting magnitude of the (1-2s)f 


component of the stator phase current. 
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The averaged data for ROTOR #2.1 is shown in Table 
4-10. Once again the (1-4s)f, (1+2s)f, and (1+4s8)f harmonic 
components were detected. Table 4-11 lists the magnitudes 
of these components for a "typical" test at each each rotor 
slip value. For all values of rotor slip, the 
non-fundamental harmonic components detected were less than 
0.25% of the fundamental component of the stator phase 
current. From Table 4-10, the variation in the magnitude of 
the (1-2s)f component is shown to be on the order of 10-25% 
of the mean value for this component. Figures 4-12 and 4-13 
show the averaged data and raw data for the (1-2s)f 
component. Figure 4-14 shows the averaged data for the 
fundamental component. Figure 4-15 shows the "typical" 
frequency spectrum for various rotor slip values. 

Tables 4-12 and 4-13 show the values of the fundamental 
and (1-2s)f components of the stator phase current measured 
for ROTOR #2 and ROTOR #2.1, and the simulation results for 
a broken rotor bar case. Comparison of the data shown in 
Table 4-12 shows that the fundamental component of the 
stator phase currents for both rotors and the simulation 
are approximately the same. For the (1-2s)f component shown 
in Table 4-13, the values for ROTOR #2.1 are significantly 
greater (by a factor of 4 to B) than those for ROTOR #2. 
This result strengthens the hypothesis that  inter-bar 
currents are present and that significant current is 


flowing through the broken rotor bar on ROTOR #2. The 





simulation results for the (1-2s)f component indicates that 
current is still flowing in the broken rotor bar on ROTOR 


#2.1. 


I(1-2s)f 
SLIP VOLTAGE STD. DEV. STD. DEV. 
V rms A rms A rms 


Table 4-10. ROTOR #2.1 averaged data. 
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Table 4-11. ROTOR 42.1 "typical" harmonic frequency data. 
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Table 4-12. Broken bar results I(60 Hz) 






I(1-2s)f AMPS rms 


Table 4-13. Broken bar results І(1-2в)Ғ 
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EXPERIMENTAL RESULTS 
STATOR PHASE I(1-2s)f CURRENT 
ROTOR #2.1 


0.020 


© 
e 
Le 
c» 


20022 


0.008 


о 
© 
© 
ea 





0.00 0.01 0.02 0.03 0.04 
ROTOR SUF 


- +- one std dev 
x mean values 


Figure 4-12. ROTOR #2.1 I(1-2s)f averaged results. 
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Figure 4-14. ROTOR £2.1 I(60 Hz) averaged results. 
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Figure 4-15. ROTOR $2.1 "typical" frequency spectrum. 
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4.5 Additional Results for a "Good" Motor (ROTOR #3) 


= MÀ MM аы Н 


In order to see the effects on the stator phase current 
due to manufacturing asymmetries, a second "as manufac- 
tured" rotor (ROTOR #3) was installed in the test motor. 
The averaged data for ROTOR #3 is shown in Table 4-14. Once 
again the (1-4s)f, (1+2s)f, and (1+48)f harmonic components 
were detected. Table 4-15 lists the magnitudes of these 
components for a "typical" test at each each rotor slip 
value. For all values of rotor slip, the non-fundamental 
harmonic components detected were less than 0.05% of the 
fundamental component of the stator phase current. From 
Table 4-14, the variation in the magnitude of the (1-2s)f 
component is shown to be on the order of 20-30% of the mean 
value for this component. Figures 4-16 and gi show the 
averaged data and raw data for the (1-2s)f component. 
Figure 4-18 shows the averaged data for the fundamental 
component. Figure 4-19 shows the "typical" frequency 
spectrum for various rotor slip values. 

Comparison of the data for ROTOR #1 (Table 4-3) and for 
ROTOR #3 (Table 4-14) shows that the fundamental component 
of the stator phase currents for both motors is 
approximately the same. For the (1-2s)f component, the 
values for ROTOR #1 are approximately two times greater 
than those for ROTOR #3. This result indicates that the 
degree of manufacturing asymmetries in ROTOR #1 is somewhat 


larger than those in ROTOR #3. 
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Table 4-14. ROTOR #3 averaged data. 
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Table 4-15. ROTOR #3 "typical" harmonic frequency data. 
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Figure 4-16. ROTOR #3 I(1-2s)f averaged results. 
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Figure 4-17. ROTOR #3 I(1-2s)f "raw" data. 
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Figure 4-18. ROTOR #3 I(60 Hz) averaged results. 
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Figure 4-19. ROTOR #3 "typical" frequency spectrum. 
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4.6 Summary 


The fundamental and (1-2s)f harmonic components of the 
stator phase current were recorded for an "off-the-shelf" 
3-HP squirrel-cage induction motor with various rotors 
installed. Two "as manufactured" rotors with no broken 
rotor bars and a third rotor with a single rotor bar 
deliberately open-circuited were used for the experiments. 
The broken bar rotor was initially tested with only one end 
of the broken bar disconnected from the end ring. 
Subsequently, the opposite end of the broken bar was 
disconnected from the other end ring and a series of 
experiments were completed using this rotor. 

For each rotor tested, in addition to the fundamental 
and (1-2s)f components of the stator phase current, 
Significant components at the frequencies (1-4s)f, (1+2s)f, 
and (1+4s)f меге detected. The (1-4s)f and (1%4-)Ғ 
components can be produced by the third and fifth 
time-harmonic components of the stator phase currents as 
shown in Appendix D. The (1+2s)f component can be 
attributed to the torque. pulsations and the resulting speed 
oscillations at  twice-the-slip frequency [7,8,9,10,11]. 
This component can also be produced by the third 
time-harmonic of the stator phase current. In general, the 
magnitudes of the (1-4s)f and (1+4s)f components were on 
the order of one-half the magnitude of the (1-2s)f 


component. The magnitude of the (1+2s)f component ів 
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approximately the same as the (1-2s)f component. For the 
two "good" rotors (ROTOR #1 and ROTOR #3) and the broken 
bar rotor with one end disconnected (ROTOR #2), the 
Magnitude of these harmonic components was less than 0.05% 
of the fundamental component. For the broken rotor bar with 
both ends disconnected (ROTOR #2.1), the magnitude of these 
components increased five-fold. 

Comparison of the experimental and simulation results 
for a motor with and without a broken rotor bar has shown 
that a more complex simulation model is required to 
accurately predict the (1-2s)f component of the stator 
phase current. In particular, the model must be able to 
account for manufacturing asymmetries which are present to 
some degree in all induction motors and it also must be 
able to model the effects of inter-bar currents. 

The existence of inter-bar currents in a squirrel-cage 
induction motor of the type tested in this thesis 
effectively "masks" the effects of a eden. rotor bar. 
These currents reduce the magnitude of the  (1-2s)f 
component of the stator phase current to a value on the 
order of those present in a "good" rotor due to 
manufacturing asymmetries. Thus, distinguishing a broken 
rotor bar from a manufacturing asymmetry using only the 
(1-2s)f component of the stator phase current appears very 
difficult. This effect is illustrated in Figure 4-20, where 


the values of the (1-2s)f component for the four rotors 
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Figure 4-20. Simulation and experimental results I(1-2s)f. 
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tested and for the simulation program are plotted as a 
function of rotor slip. A distinction between a broken 
rotor bar and a manufacturing asymmetry is only possible 
after both ends of the broken rotor bar are disconnected 
from the end ring. 

As a final note, if the single-test diagnostic 
algorithm of Kliman et al. [7] (see Chapter 1) were used, a 
motor with either ROTOR #1, ROTOR #2, or ROTOR #3 installed 
would be declared a "good" motor. For a motor with ROTOR 
#2.1 installed, the measured dB value (-54 dB) falls in the 
region between a "good" motor and a "broken rotor bar" 


motor and thus would be declared "suspect". 
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CHAPTER 5 


CONCLUSIONS AND RECOMMENDATIONS FOR FUTURE RESEARCH 


5.1 Summary of Thesis 


The purpose of this thesis was to critically evaluate 
one method for the detection of broken rotor bars in 
induction motors using stator current and voltage 
measurements. This effort supported the ongoing development 
of a failure analysis system for electric machines by the 
Ии: Laboratory for Electromagnetic and Electronic 
Systems. The hypothesis of the method, given a sinusoidally 
applied voltage, was that the presence of certain harmonics 
in the stator currents could be used to detect the presence 
of broken rotor bars. 

To support the evaluation, a system of first-order 
differential equations describing the electrical 
performance of a three-phase induction motor was developed 
using stator phase currents and rotor loop currents as 
state variables. Each stator phase and rotor loop was 
described in terms of a resistance and inductance. The 
stator phase voltages (assumed to be sinusoidal) were the 
driving functions for the systen. 

A relationship between the standard single-phase, 
equivalent-circuit model for an induction motor and the 
system of equations was derived (Equations 2-100 to 2-106). 
With this relationship, the electrical parameters needed to 


solve the system of equations were easily calculated from 
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the single-phase, equivalent-circuit parameter values which 
are often available from the manufacturer, or are easily 
measured in the laboratory. This relationship was useful in 
that the solution to the system of equations could be 
obtained without disassembling the motor, removing the 
rotor, and measuring the geometrical and material 
properties of the motor such as the rotor radius, the 
alrgap length, and conductivity of the rotor bars. Thus, by 
performing the standard no-load, locked-rotor, and DC tests 
and knowing the number of rotor bars, the system of 
equations could be used to simulate a three-phase induction 
motor with an arbitrary number of rotor bars. 

A FORTRAN simulation routine was developed to solve the 
system of equations describing the electrical performance 
of an induction motor with or without a broken rotor bar. 
The routine uses a fixed-step, fourth-order Runge-Kutta 
procedure to numerically integrate the system of equations. 
A matrix inversion routine, Cholesky’s method with partial 
pivoting, was required to invert the time-varying 
inductance matrix for each value of time. In addition to 
the FORTRAN simulation program, two PRO-MATLAB routines 
were also employed. The PRO-MATLAB eigenvalue routine was 
used to aid in determining the fixed time step for ae 
Runge-Kutta integration. A PRO-MATLAB FFT routine was used 


to transform the output data of the FORTRAN program from 
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the time domain to the frequency domain. Thus the 
fundamental and (1-2s)f frequency components of the stator 
phase currents could be determined. 

The FORTRAN simulation program and PRO-MATLAB FFT 
routine were initially validated using a hypothetical 3-bar 
rotor induction motor. The errors between the simulation 
results and exact solutions for both the no-broken-bar and 
the one-broken-bar cases were less than 1%. The simulation 
program accurately solved the system of equations derived 
for the stator phase currents and rotor loop currents. 
Thus, the simulation permitted numerical experiments which 
were used to test the failure detection hypothesis. 

For the numerical experiments, the single-phase, 
equivalent-circuit parameter values for an "off-the-shelf" 
3-HP squirrel-cage induction motor were used as inputs to 
the simulation routine. Thus, the simulations represented 
the operation of the same motors which were used for the 
physical experiments. The results of the simulation for a 
no-broken-bar case showed only a fundamental component of 
the stator phase currents. However, for the one-broken-bar 
case, the simulation showed that an additional component at 
a frequency of (1-2s)f was also present in the stator phase 
current. The magnitude of this component was approximately 
2% of the magnitude of the fundamental component for the 
full load operating speed of the motor. Thus, the failure 


detection hypothesis appeared plausible and physical 
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experiments were undertaken to verify that the (1-2s)f 
component of the stator phase current could be used to 
detect the presence of a broken rotor bar. 

The physical experiments were conducted with 
commercially-available induction motors. The fundamental 
and (1-2s)f harmonic components of the stator phase current 
were recorded for an "off-the-shelf" 3-HP squirrel-cage 
induction motor with various rotors installed. Two "as 
manufactured” rotors with no broken rotor bars and a third 
rotor with a single rotor bar deliberately disconnected 
from the end ring by removing a small section of the rotor 
bar were used for the experiments. The broken-bar rotor was 
initially tested with only one end of the broken bar 
disconnected from the end ring. Subsequently, the opposite 
end of the broken bar was cut away from the other end ring 
and a series of experiments was completed using this rotor. 

For the two “good” rotors tested, the (1-2s)f component 
of the stator phase current was dosssed. The magnitude of 
this component was approximately 0.05% of the fundamental 
component. The (1-2s)f component was present in a "good" 
rotor due to tan asymmetries which are present to 
some degree in all machines. The measured values of the 
(1-2s)f component for the broken bar rotor with one end 
disconnected were also only 0.05% of the fundamental 


component. Thus, using only the (1-2s)f component of the 


stator phase current, the distinction between manufacturing 
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asymmetries and a broken rotor bar was impossible. This 
result was significantly different than that predicted by 
the simulation program. 

In order to explain the differences between the 
simulation and experimental results, a hypothesis that 
current was still flowing through the broken rotor bar was 
investigated. One possible mechanism was that the bar was 
not completely broken. This could be the case in the 
present experiments if the rotor bar was not completely 
disconnected from the end ring during milling so that path 
for current to flow still exists. After removing the rotor 
and visually inspecting the broken rotor bar, it was 
concluded that the rotor bar was completely cut away from 
the end ring. Thus it appeared certain that current is not 
flowing into the broken bar across the Creak fron the end 
ring. 

Another mechanism that could result in a current 
flowing through a broken rotor bar was that of inter-bar 
currents. Inter-bar currents are currents which flow 
between a rotor bar and adjacent rotor bars through the 
rotor laminations. The presence of such currents could be 
expected to "short circuit” the effect of a break in a 
rotor bar, hence resulting in a current flow in the broken 
bar that would not be predicted by the simulation program. 

In order to investigate the hypothesis that inter-bar 


currents were masking the effects of the break in the rotor 
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bar during experiments with the broken bar rotor, the 
opposite end of the broken rotor bar was also cut away from 
the end ring. This was expected to result in a further 
decrease in the magnitude of the current flowing through 
the broken rotor Баг, thus increasing the resulting 
magnitude of the (1-2s)f component of the stator phase 
current. The measured values of the (1-2s)f component for 
this rotor were on the order of 4 to 8 times greater than 
the previous results. This strengthened the hypothesis that 
inter-bar currents were flowing and that significant 
current was flowing through the broken rotor bar with one 
end disconnected during the previous tests. In addition, 
the simulation results were still an order of magnitude 
greater than the measured results for the broken rotor bar 
with both ends disconnected, thus indicating that current 


was still flowing in the broken rotor bar. 


5.2 Conclusion 


The existence of inter-bar currents in a squirrel-cage 
induction motor of the type tested in this thesis 
effectively "masks" the effects of a broken rotor bar. 
These currents reduce the magnitude of the (1-2s)f 
component of the stator phase current to a value on the 
order of those already present in a "good" rotor due to 
manufacturing asymmetries. The distinction between a broken 
rotor bar and a manufacturing asymmetry is only possible 


after both ends of a rotor bar are cut away from the end 
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ring. Thus, distinguishing a broken rotor bar from a 
manufacturing asymmetry using only the (1-2s)f component of 


the stator phase current appears highly improbable. 


5.3 Recommendations for Future Work 


Although the results of this thesis show that the 
detection of a broken rotor bar in a squirrel-cage 
induction motor of the type tested is highly improbable 
using only the (1-2s)f component of the stator’ phase 
current, examination of the higher-order stator phase 
current harmonics may provide a basis for distinguishing 
between a broken rotor bar and manufacturing asymmetries. 
Kliman et al. [7] has shown that the higher-order harmonic 
components of the airgap flux resulting from a broken rotor 
bar are typically much larger than those resulting from 
manufacturing asymmetries. Thus, an investigation of the 
effects on the higher-order stator phase current harmonics 
due to a broken rotor bar is needed if the failure analysis 
system is to be based solely on the measurements of 
terminal voltages and currents. In addition, the effects of 
inter-bar currents would have to be modeled in detail. 

In order to use the simulation program developed in 
this thesis in a failure analysis system for electrical 
machines, the program must be modified to include the 
effects of inter-bar currents. This requires an in-depth 
analysis on the existence of and effects of inter-bar 


currents in squirrel-cage induction motors. In addition to 
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modifying the simulation program, changes in the current 
manufacturing process for squirrel-cage rotors may be 
recommended based on this analysis (e.g., insulation 
between the rotor bars and laminations). 

Axial flux monitoring [11] is another method which has 


been successful in detecting broken rotor bars, as well as 


various other faults, in operating induction machines. 
Although this method requires that an external sensing coil 
be attached to a machine, it appears to be very promising 
and can be implemented into a failure analysis system. This 
method is currently being investigated at the M.1.T. 


Laboratory for Electromagnetic and Electronic Systems [21]. 
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APPENDIX A 


ELEMENTS FOR SYSTEM MATRICES 


The following sections give the expressions for the 


elements of the voltage, inductance, and effective 
resistance matrices derived in Chapter 2 (NOTE: m=row, 
n=column). Refer to the List of Symbols given in section 2.2 


for a definition of symbols. 
A-1 Voltage Matrix [v] 


for Wel sags 





YS (дуе Л) 
vav, (АС?) 
Va = V. CES) 
for mz4,5,...N.  +3 
v. = 0 f СЕ 
A-2 Inductance Matrix [L] 
Xue HBSEHYS 55 
Bl Zn (A) 
where l, is given by equation 2-26. 
for men m,n=1,2,3 
me (A-6) 
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forn mel, 2,3 = 4 5, М3 


LL, Main Ын. w,t+0,)-(m-1)22) CAT) 
\ RB Э 


where M is given by equation 2-52. 


fiorngmzi,55..N»t3 nzl,2;29 


2(m-3)-1)r 
Lea Msin{ p( CONSE, uc aeg) qu i27 (AR) 
floraMm-nz41,5;-..-Nimnt393 
PUB vs -5;* Lina (99 


a ctc егам У а Е) 


where [, 18 given by equation 2-42. 


for nzm+1 m=1,5,...N,. F +3 


-L; 


last.) unas (A-10) 
Move: for m=N\,,+3, n=4 В 1910) 
Бот al Am Ed O. NS 
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Ls cr m- (A ИЮ 
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for 


for 


nsm,m+1l,m-] п,шй=4‚,5,...Мь+3 
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шо me) A) 
Effective Resistance Matrix [RR] 

m=n=1,2,3 

КЕ и =R, CASES) 
men m,n=1,2,3 

КЕ „ „20 (А - 14) 
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ОК 
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(note: this expression assumes the mechanical speed of 
the rotor is a constant) 


BERN UN n n=162.3 


mI 


RR,,7p«o,Mcos| p 


wt+0,)-(n-1) 2) (A-16) 


(note: this expression assumes the mechanical speed of 
the rotor is a constant) 
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ЕН = 0 (A-20) 
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APPENDIX B 


SIMULATION PROGRAMS AND DATA FILES 


B-1 FORTRAN SIMULATION PROGRAM 


The FORTRAN simulation program can be used to solve the 
system of equations describing the electrical performance of 
an induction motor derived in Chapter 2 for the time-varying 


Stator phase and rotor loop currents. A  fixed-step, 
fourth-order Runge-Kutta procedure is used to numerically 
integrate the system of equations. The time-varying 


inductance matrix is inverted for each value of time using 
Cholesky’s method with partial pivoting. The program can be 
used to simulate an induction motor with either no broken 
rotor bars or one broken rotor bar. A flowchart for the 
program is shown in Figure 3-1. The program is divided into 
six parts, a main program and five’ subroutines. The 
following paragraphs describe the purpose of each part of 
the simulation program. 


Main program: The main program performs several 
functions. The primary function is to act as a buffer 
and pass data between the subroutines. In addition, the 
main program keeps track of the simulation time, 
calculates the value of the current vector at the end 
of each time step, and writes the results to an 
external file. 


Subroutine INPUT: This subroutine is called once by the 
main program. It performs two functions. First, it 
reads the data contained in an external input file. The 
external input file provides a description of the 
machine being simulated. It includes the single-phase, 
equivalent-circuit model values, number of rotor bars, 
and number of pole-pairs for the machine, the rotor 
slip, and the supply voltage and frequency. The 
Simulation time parameters ( i.e., step-size and stop 
time) and a broken bar flag are also included in the 
input file. Second, using the relations derived in 
section 2.3, the coefficients for the inductance and 
effective resistance matrix elements are calculated and 
passed to the main program. 


Subroutine MOTOR: This subroutine is called for each 
value of time. It performs two functions; it calculates 
the elements of the voltage, inductance, and effective 
resistance matrices, and it reduces the system of 
equations. For the voltage vector, a balanced 
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three-phase supply is assumed. The inductance and 
effective resistance elements are calculated using the 
relations given in Appendix A. For both a “no broken 
bar" simulation and a "one broken bar" simulation, the 
system of equations can be reduced from N,4*3 to М,,%2 
since the sum of the rotor loop currents equals zero. 
This condition is imposed on the system of equations by 
eliminating rotor loop Na current. Column Na+3 із 
subtracted from columns 4 through Nat2 and column Nast3 
is deleted from the inductance and effective resistance 
matrices. Next, row N,a+3 is deleted from each matrix 
to complete the reduction. In addition, for a "one 
broken bar" simulation (assume rotor bar Мљ-1 is 
broken), the system of equations can be reduced again 
by setting rotor loop N&-1 current equal to rotor loop 
Na-2 current. For this condition, column N,,+2 is added 
to column Naat] and column N,,*2 is deleted from the 
inductance and effective resistance matrices. Faraday’s 
law is satisfied by adding row N&4*2 to row N,*1 and 
deleting row N4*2 from each matrix to complete the 
reduction. 


Subroutine CURRENT: This Bubroutine is used to 
calculate the intermediate values of the current vector 
for the Runge-Kutta method. 


[ESTI- 11 O 


where: 


([EST]= intermediate value of current vector 


[1]= initial current vector 


Subroutine DRIVE: The purpose of this subroutine is to 
calculate the "driving function" for each value of time 
and current, i.e., 


driving function- [V]- [RR )(1] 


Subroutine CHLSKY: This subroutine performs the 
inductance matrix inversion to calculate the derivative 
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of the current vector for each time value. For the 
fourth-order Runge-Kutta procedure, this calculation is 
performed four times per time step. 
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PROGRAM LISTING 


PROGRAM INDUCTION MOTOR 


FILENAME: MARK.FOR 
AUTHOR: MARK S. WELSH 
REFERENCES: CHAPTERS 2 AND 3 OF THESIS 
SECTIONS 3-4 AND 6-5 OF REFERENCE [13] 


THIS PROGRAM SIMULATES THE ELECTRICAL OPERATION OF A 
3-PHASE SQUIRREL-CAGE INDUCTION MOTOR WITH NRB ROTOR 
BARS. THE PROGRAM SOLVES THE SYSTEM OF EQUATIONS IN 
THE FORM OF 

[V]=[L]d/dt[1]+[RR](1] 
FOR THE STATOR PHASE AND ROTOR LOOP CURRENTS. A 
FIXED-STEP, FOURTH-ORDER, RUNGE-KUTTA INTEGRATION 
ROUTINE IS USED. IN ADDITION CHOLESKY’S METHOD WITH 
PARTIAL PIVOTING IS USED TO SOLVE FOR THE VALUES OF 
d/dt[1] FOR EACH TIME VALUE (I.E., INVERTS [L) FOR 
EACH TIME VALUE). THE PROGRAM IS BROKEN DOWN INTO 
SIX SECTIONS (MAIN PROGRAM AND FIVE SUBROUTINES). 
THE PROGRAM WILL ALLOW FOR SIMULATION OF EITHER NO 
BROKEN ROTOR BARS OR ONE BROKEN ROTOR BAR. 


LIST OF VARIABLES FOR MAIN PROGRAM 


NAMOUT OUTPUT FILE NAME 

ВКХ FLAG FOR A BROKEN ROTOR BAR 

NRB NO. OF ROTOR BARS 

k NO. OF SYSTEM EQUATIONS (NRB+3) 

Kl NO. OF SYSTEM EQUATIONS FOR THE NO BROKEN 
BAR CASE (NRB+2) 

Por NO. OF SYSTEM EQUATIONS FOR THE ONE BROKEN 
BAR CASE (NRB+1) 

LS STATOR PHASE SELF INDUCTANCE 

LSL STATOR PHASE LEAKAGE INDUCTANCE 

RS STATOR PHASE RESISTANCE 

MSR STATOR PHASE-ROTOR LOOP MUTUAL INDUCTANCE 
COEFFICIENT 

LR ROTOR LOOP SELF INDUCTANCE 

LRLB ROTOR BAR LEAKAGE INDUCTANCE 

RRB ROTOR BAR RESISTANCE 

E INPUT VOLTAGE AMPLITUDE 

W INPUT VOLTAGE ELECTRICAL FREQUENCY 

WM ROTOR MECHANICAL SPEED 

T SIMULATION TIME 

INC TIME STEP FOR INTEGRATION 

TSTOP SIMULATION STOP TIME 

TOUT SIMULATION TIME TO BEGIN WRITING OUTPUT 

FACTOR INTERMEDIATE TIME STEP FOR INTEGRATION 

TIME INTERMEDIATE SIMULATION TIME 

EP NO. OF MACHINE POLE-PAIRS 

AMPS CURRENT VECTOR 
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ACERCAS SEE) 


() 2) ооо 


QOO 


соо 


10 


ооо 


v VOLTAGE VECTOR 


L INDUCTANCE MATRIX 

R EFFECTIVE RESISTANCE MATRIX 

VRI DRIVE VECTOR (VRIz[V]-[RR][I]) 

EST INTERMEDIATE CURRENT ESTIMATE 

Rk 1 INITIAL VALUE OF d/dt[I] FOR EACH STEP 
RK2 INTERMEDIATE VALUE d/dt[I} FOR EACH STEP 
ЕКЗ INTERMEDIATE VALUE d/dt[I] FOR EACH STEP 
RK4 LAST VALUE OF d/dt[1} FOR EACH STEP 


VARIABLE DECLARATIONS 


CHARACTER*20 NAMOUT 

INTEGER BKN,NRB,K,K1,K2 

DOUBLE PRECISION LS,LSL,RS,MSR,LR,LRLB,RRB,E,W,WM,T, 
* INC, TSTOP, TOUT, FACTOR, TIME, PP, AMPS(100,1),V(100,100), 
*L(100,100),R(100,100),VRI(100,1),EST(100,1), 
*RK1(100,1),RK2(100,1),RK3(100,1),RK4(100,1) 


GET INPUT PARAMETERS 


CALL INPUT(LS,LSL,RS,MSR,LR,LRLB,RRB,BKN,E,W,WM,NRB, 
*TSTOP,TOUT,INC,PP) 


SET UP OUTPUT FILE 


PRINT *,'ENTER OUTPUT FILE NAME’ 
READ *,NAMOUT 
OPEN(UNITz7,FILEzNAMOUT,STATUS-z 'UNKNOMN') 


DETERMINE NUMBER OF SYSTEM EQUATIONS 


K=NRB+3 
КІ-К-1 
IF(BRN.LE.O) THEN 
К2-К1 
ELSE 
K2=K1-1 
ENDIF 


INITIALIZE SIMULATION TIME AND CURRENTS 
РО о = оке 
AMPS(J,1)=0. 
CONTINUE 
T20. 
SET UP SYSTEM MATRICES FOR TIME=T 


CALL MOTOR(V,L,R,K,K1,K2,T,LS,LSL,RS,MSR,LR,LRLB,RRB, 
*BKN,E,W,WM,NRB, PP) 


START SIMULATION LOOP 
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WRITE OUTPUT TO DATA FILE (TIME, PHASE A CURRENT, 
AND ROTOR LOOP 1 CURRENT) 


ТЕСТ СЕ ТОСТ) THEN 
WRITE(7,999)T,AMPS(1,1),AMPS(4,1) 

ENDIF 

CALCULATE RK1 (d/dt[1]) FOR TIME=T 

CALCULATE VRI FOR TIME=T, I=AMPS 

CALL DRIVE(R,K2,AMPS,V,VRI) 

CALL CHLSKY(L,K2,K2+1,VRI,RK1) 

CALCULATE RK2 (d/dt[{1])) FOR TIME=T+INC/2 


TIME=T+INC/2.0 
FACTOR=INC/2.0 


ESTIMATE CURRENT FOR TIME=T+INC/2 USING RK1 
CALL CURRENT(AMPS,RK1,FACTOR,K2,EST) 
SET UP SYSTEM MATRICES FOR TIME=T+1NC/2 


CALL MOTOR(V,L,R,K,K1,K2,TIME,LS,LSL,RS,MSR,LR,LRLB, 
*RRB,BKN,E,WM,NRB, PP) 


CALCULATE VRI FOR TIME=T+INC/2, I=EST 
CALL DRIVE(R,K2,EST,V,VRI) 

CALL CHLSKY(L,K2,K2+1,VRI,RKĶK2) 

CALCULATE RK3 (d/dt[I]) FOR TIME=T+INC/2 
ESTIMATE CURRENT FOR TIME=T+INC/2 USING RK2 
CALL CURRENT(AMPS,RK2,FACTOR,K2,EST) 
CALCULATE VRI FOR TIME=T+INC/2, I=EST 
CALL DRIVE(R,K2,EST,V,VRI) 

CALL CHLSKY(L,K2,K2+1,VRI,RK3) 

CALCULATE RK4 (d/dt[{1]) FOR TIME=T+INC 
TIME=T+INC 


ESTIMATE CURRENT FOR TIME=T+INC USING RK3 
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CALL CURRENT(AMPS,RR3,INC,K2,EST) 
SET UP SYSTEM MATRICES FOR TIME=T+INC 


CALL MOTOR(V,L,R,K,K1,K2,TIME,LS,LSL,RS,MSR,LR,LRLB, 
*RRB,BKN,E,W,WM,NRB, PP) 


CALCULATE VRI FOR TIME=T+INC, I=SEST 

CALL DRIVE(R,K2,EST,V,VRI) 

CALL CHLSKY(L,K2,K2+1,VRI,RK4) 

CALCULATE NEW CURRENT VALUES AND INCREMENT TIME 

DO 20 J=1,K2 
AMPS(J,1)=AMPS(J,1)+(INC/6.0)*(RK1(J,1)+2#RK2(J,1)+ 

t 2*RK3(J,1)+RK4(J,1)) 

CONTINUE 

T=TIME 

CONTINUE SIMULATION IF T<TSTOP ELSE STOP 

IF(T.LE.TSTOP) GO TO 1 


CLOSE(UNIT=7) 
STOP 


FORMAT(3(3X,F12.7)) 
END 
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SUBROUTINE INPUT(LS,LSL,RS,MSR,LR,LRLB,RRB,BKN,E,W, 
*WM,NRB, TSTOP, TOUT, INC, PP) 


THIS SUBROUTINE IS THE INPUT INTERFACE FOR THE 
SIMULATION PROGRAM. REQUIRED INPUT VALUES ARE READ 
IN FROM AN EXISTING INPUT DATA FILE AND THE VALUES 
OF INDUCTANCES,RESISTANCES,ETC. ARE CALCULATED FROM 
THE RELATIONS DERIVED IN SECTION 2.3 AND PASSED TO 
THE MAIN PROGRAM 


LIST OF VARIABLES 


NAMEIN INPUT DATA FILENAME 
STEPS NO. OF TIME STEPS PER SECOND 


5 ROTOR SLIP 

PI RADIANS IN A HALF CIRCLE 

h CONSTANT FACTOR FOR CALCULATIONS 

R1 CIRCUIT MODEL STATOR PHASE RESISTANCE 

L1 CIRCUIT MODEL STATOR PHASE LEAKAGE INDUCTANCE 
L12 CIRCUIT MODEL MUTUAL INDUCTANCE 

L2 CIRCUIT MODEL ROTOR LEAKAGE INDUCTANCE 

R2 CIRCUIT MODEL ROTOR RESISTANCE 


VARIABLE DECLARATIONS 


CHARACTER *20 NAMEIN 

INTEGER BKN,NRB,PP,STEPS 

DOUBLE PRECISION TSTOP,TOUT,LS,LSL,RS, MSR, LR, 
*LRLB,RRB,E,W,WM,S,PI,INC,K,R1,L1,L12,L2,R2 


READ IN INPUT DATA FILE 


PRINT *,'ENTER INPUT FILENAME' 

READ *,NAMEIN 
OPEN(UNIT=8,FILE=NAMEIN, STATUS=’ UNKNOWN’ ) 
READ(8,*)R1,L1,L12,L2,RZ 

READ(8,* )NRB, BEN 

READ(8,*)E,W,S,PP 
READ(8,*)STEPS,TSTOP, TOUT 

CLOSE(UNIT=8) 


DEFINE PARAMETERS 


PI=4*ATAN(1.0) 

K=PI*+PP/NRB 

LSL=L1 

RS=Rl 

ава 12730 

MSR:8*SIN(K)*L12/(3*PI) 
LRLB=(L2-(K**2/(SIN(K)*¥%2)-1)*L12)%4*NRB/(3*PI1**#2) 
RRB=4*NRB/(3*P1*4*2)*R2 

LR=16%(NRB-1) *PP*#2%*L12/(3*NRB#* #2) 
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PNC=1.0/STEPS 
WM=Wx (1-S)/PP 
RETURN 

END 


АЫ 
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SUBROUTINE MOTOR(V2,L2,R2,K,K1,K2,TIME,LS,LSL,RS,MSR, 
*LR,LRLB,RRB,BKN,E,W,WM,NRB, PP) 


THIS SUBROUTINE CALCULATES THE ELEMENTS IN 

EACH MATRIX OF THE SYSTEM OF EQUATIONS FOR A GIVEN 
VALUE OF TIME. THE MATRICES ARE SET-UP AND THEN 
REDUCED FOR THE CONDITION THAT THE SUM OF THE ROTOR 
LOOP CURRENTS IS ZERO. FOR A BROKEN BAR THE SYSTEM IS 
AGAIN REDUCED BY ONE SINCE THE ROTOR LOOP CURRENTS 
ASSOCIATED WITH A BROKEN BAR ARE SET EQUAL. FOR THIS 
SIMULATION ROTOR BAR NRB-1 IS THE BROKEN BAR. 


LIST OF VARIABLES 

V1,R1,L1 MATRICES V,R,L REDUCED FOR SUM OF ROTOR 
LOOP CURRENT=0 CONDITION 

V2,R2,L2 MATRICES V1,R1,L1 REDUCED FOR BROKEN ROTOR 
BAR CONDITION 

VARIABLE DECLARATIONS 


INTEGER k,K1,kK2,BKN,NRB, PP 
DOUBLE PRECISION V(100,1),R(100,100),L(100,100), 


*V1(100,1),R1(100,100),L1(100,100),V2(100,1), 


28201005100) ,22(100,100), TIME,LS,LSL,RS,MSR,LR,LRLB, 


*RRB,E,W,WM,PI 


PI=4#ATAN(1.0) 
CALCULATE VOLTAGE VECTOR 


DO 10 J=1,k 
IF(J.GT.3) THEN 


V(J,1)=0.0 
ELSE 
V(J,1)=E*COS(W*TIME-(J-1)*2*P1/3.0) 
ENDIF 
CONTINUE 


CALCULATE INDUCTANCE AND EFFECTIVE RESISTANCE MATRIX 
ELEMENTS 


DO ZO 1-1-К 
DO 30 J=zI,K 
IF(I.LE.3.AND.J.LE.3) THEN 
IF(I.EQ.J) THEN 
L(I,J)=LS+LSL 


R(I,J)=RS 
ELSE 
L(I,J)z-LS/2.0 
R(I,J)=0.0 
ENDIF 


180 





30 
20 


IIS 


OD» 


aca 


60 


ELSENHUT.LE.3 AND. J.CT.3) THEN 
L(1,J)=MSR*SIN( PP*WM*TIME+(2*(J-3)-1)# 
PI/NRB-(1I-1)*21P1/3.0) 
R(I,J)zPP*WM£MSRECOS(PP£WMETIME*(2t(J-3) 
-1)ЖРІ/МЕВ-(1-1)%2%Р1/3.0) 
ELSE 
IF(I.EQ.J) THEN 
L(I,J)-LR*2*LRLB 
R(I,J)z2*RRB 
ELSEIF(J.EQ.(1+1))THEN 
L(I,J)z-LR/(NRB-1)-LRLB 
R(I,J)=-RRB 


ELSE 
L(I,J)=-LR/(NRB-1) 
R(I,J)=0.0 

ENDIF 

ENDIF 


R(J,I)=R(I,J) 
L(J,I)=L(I,J) 
CONTINUE 
CONTINUE 
L(4,K)=-LR/(NRB-1)-LRLB 
L(K,4)=L(4,k) 
R(4,K)=-RRB 
R(k,4)=R(4,K) 


REDUCE EQUATIONS FOR SUM OF ROTOR LOOP CURRENTS 
EQUAL ZERO CONDITION 


DO 40 I=1,k1 
DO 50 J=1,F1 
ШЕИ ПА ЕУ НЕМ 
R1(1,J)=R(1,J) 
L1(1,J)=L(1,J) 
ELSE 
R1(I,J)=R(I,J)-R(I,K) 
ШЕГШЕ ЕЕ ШЕЕ = ШЫК) 
ENDIF 
CONTINUE 
V1(I,1)zV(I,1) 
CONTINUE 


FOR A BROKEN BAR REDUCE EQUATIONS AGAIN ELSE 
SET V2=V1,R2=R1,L2=L1 AND RETURN TO MAIN PROGRAM 


IF(BKN.LE.O) GO TO 1000 

REDUCING EQUATIONS FOR BROKEN BAR 

DO 60 I=z1,K1 
R1(I,K2)=R1(I,K2)+R1(I,K1) 


ПИКЕ ЕТТЕ) ТКА) 
CONTINUE 
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DO 70 I-1,K1 
R1(K2,1)=R1(K2,1)+R1(K1,1) 
L1(K2,1)=L1(K2,1)+L1(K1,1) 

CONTINUE 


SETTING FINAL VALUES FOR V,R,L MATRICES 


DO 80 I=1,K2 
DO 90 J=1,K2 
R2(I,J)=R1(I,J) 
EZI E Л) 
CONTINUE 
CONTINUE 
RETURN 
END 
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SUBROUTINE CURRENT(A,B,C,N,RESULT) 


C 
C THIS SUBROUTINE CALCULATES THE INTERMEDIATE 
€ VALUES OF CURRENT FOR THE INTEGRATION ROUTINE. 
С ([I(EST) J=[I(start) ]+d/dt[I]tdt 
E 
C LIST OF VARIABLES 
€ 
C N NO. OF EQUATIONS 
C A INITIAL VALUE OF CURRENTS 
C B d/dt[I] 
€ E dt (INCREMENTAL TIME) 
Ç RESULT ESTIMATED CURRENT (A«*B*C) 
C 
C VARIABLE DECLARATIONS 
С 

INTEGER N 

DOUBLF PRECISION A(100,1),B(100,1),RESULT(100,1),C 
e 
C CALCULATE CURRENT ESTIMATE 
E 

BO 10. J=1,N 

RESULT(J,1)=A(J,1)+B(J,1)*C 

10 CONTINUE 

RETURN 

END 
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SUBROUTINE DRIVE(A,K2,B,C,RESULT) 


THIS SUBROUTINE CALCULATES THE VALUE OF 
[V]-[RR])[1] FOR EACH VALUE OF TIME. 


LIST OF VARIABLES 


A EFFECTIVE RESISTANCE MATRIX [RR] 
B CURRENT VECTOR [T] 
E VOLTAGE VECTOR [V) 


RESULT — [V]-([RR][1] 
VARIABLE DECLARATIONS 


INTEGER K2 
DOUBLE PRECISION A(100,100),B(100,1),C(100,1), 


*RESULT(100,1) 


CALCULATE RESULT 


DO 10 I=1,K2 
RESUÚLT(I,1)=0.0 


DETERMINE [RR][I] 
DO 20 J=1,K2 
RESULT(1,1)=RESULT(1,1)+A(1,J)*B(J,1) 
CONTINUE 
DETERMINE [V]-[RR][1] 
RESULT(1,1)=C@l,1)-WESUET(1,1) 
CONTINUE 


RETURN 
END 
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SUBROUTINE CHLSKY(A1,N,M,B,RESULT) 


THIS SUBROUTINE CALCULATES THE VALUES OF d/dt[I] 

FOR EACH VALUE OF TIME USING CHOLESKY’S METHOD WITH 

PARTIAL PIVOTING. THE EQUATIONS ARE IN THE FORM: 
[L]d/dt[I]z[V)-[RR][I] 


LIST OF VARIABLES 


A WORKING MATRIX [L;V-RR I] 
A1 INDUCTANCE MATRIX [L) 
B DRIVE VECTOR [V]-I[RRJII] 


RESULT CURRENT DERIVATIVES d/dt[I] 

PIVOT VALUE OF PIVOT ELEMENT 

SWAP DUMMY VARIABLE FOR CHANGING ROWS 

SUM VARIABLE FOR MULTIPLICATION OPERATIONS 


VARIABLE DECLARATIONS 


INTEGER N,M 
DOUBLE PRECISION A(100,101),A1(100,100),B(100,1), 
*RESULT(100,1),PIVOT,SWAP,SUM 


SET UP A MATRIX [A]z[A1;|B] 
(ADD DRIVE COLUMN TO INDUCTANCE MATRIX) 


DO 1 1=1,N 
рО 2 4-1,М 
IF(J.NE.M) THEN 
Ат ОСА С) 
ELSE 
А(Т,Ј)=В(1,1) i 
ENDIF 
CONTINUE 
CONTINUE 


PERFORM PARTIAL PIVOTING 
DO 10 I=1,N 
FIND LARGEST PIVOT 


PIVOT=A[1,1) 

TAL 

DO 20 J=I+1,N 
IF(ABS(A(J,I)).LE.ABS(PIVOT)) GO TO 20 
PIVOT=A(J,1) 
IBzJ 

CONTINUE 

IF(IL.EQ.I) GO TO 10 


INTERCHANGE ROWS TO PUT PIVOT ON DIAGONAL 
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DO 30 K=1,M 
SWAP-A(I,K) 
A(I,K)=A(IL,K) 
A(IL,K)=SWAP 

CONTINUE 

CONTINUE 


CALCULATE FIRST ROW OF UPPER TRIANGULAR MATRIX 
DO 40 J=2,M 

A(1,J)=A(1,J)/A(1,1) 
CONTINUE 


DETERMINE ALL OTHER ELEMENTS OF UPPER AND LOWER 
TRIANGULAR MATRICES 


ШО 50-225 
LOWER TRIANGULAR MATRIX ELEMENTS 
DO 60 J=I,N 
SUM=0.0 


ПОШ ПОК = ЫЕ 
SUM=SUM+A(J,K)*A(K,1) 


CONTINUE 
A(J,1)=A(J,1)-SUM 
CONTINUE 


UPPER TRIANGULAR MATRIX ELEMENTS 


DO 80 J=I+1,M 
SUM=0.0 
DO 90 K=1,1-1 

SUM=SUM+A(1,K)*A(K,J) 

CONTINUE 
A(I,J)=(A(I,J)-SUM)/A(I,]I) 

CONTINUE 

CONTINUE 


USE BACK SUBSTITUTION TO FIND d/dt[TI] 


RESULT(N,1)=A(N,M) 
DO 100 I=1,N-1 
J=N-I 
SUM=0.0 
DO 110 K=J+1,N 
SUM=SUM+A(J,K) *RESULT(K,1) 
CONTINUE 
RESULT(J,1)=A(J,M)-SUM 
CONTINUE 
RETURN 
END 
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b. SAMPLE INPUT FI 





The following input file is for a 45-bar rotor with a 
broken rotor bar at s=0.04. The values are in the following 
format (i.e., refer to subroutine input above): 


equivalent-circuit values (R1,L1,L12,L2,R2) 

no. of rotor bars, broken bar flag 

input voltage, electrical frequency, slip, pole-pairs 
no. of steps/sec, stop time, start print time 


.859,.0046,.0704,.0046,.5612 
45,1 

W59.71,3/7,.04,2 
1536,4.3334,1.0 


Table B-1. Sample simulation input file. 


187 





с. SAMPLE OUTPUT FILE 


The following output file is for a 45-bar rotor with a 
broken rotor bar. This is a portion of the output file 
generated from the simulation program using the input file 
given in section 1.1.1. Plots of stator phase a current and 
rotor loop 1 current are included. 





TIME (sec) STATOR PHASE A ROTOR LOOP 1 
CURRENT (amps) CURRENT (amps) 

1.0000 9.388 2.926 
1.0007 10.963 2.902 
1.0013 11.878 2.878 
1.0020 12.080 2.854 
1.0026 1010257516 2.830 
1.0033 10.340 2.805 
1.0039 8.503 2.780 
1.0046 6.158 2.754 
1.0052 3.444 2.729 
1.0059 0,528 2.703 
1.0065 -2.424 2.676 
1.0072 -5.225 2.650 
1.0078 -7.711 2.623 
1.0085 -9.733 2.596 
1.0091 -11.169 2.569 
1.0098 -11.933 2.542 
1.0104 -11.980 2514 
ПОТОМ -11.308 2.486 
1.0117 -9.957 2.458 
1.0124 -8.008 2.430 
1.0130 -5.580 2.401 
1.0137 -2.817 2312 
1.0143 0.114 2.343 
1.0150 3.036 | 2.314 
1.0156 Б. 75 2.284 
1.0163 8.165 2.254 
1.0169 10.064 2.224 
1.0176 11959857 2.194 
1.0182 11.968 2.164 
1.0189 11.858 2.133 
1.0195 11.037 2.102 
1.0202 97592 2001 
1.0208 7.494 2.040 
1.0215 4.986 2.009 
1.0221 2.179 1.977 
1.0228 -0.758 1.945 
1.0234 -3.648 1.913 
1.0241 -6.319 1.881 


Table B-2. Sample simulation output file. 
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TIME (sec) STATOR PHASE A ROTOR LOOP 1 


CURRENT (amps) CURRENT (amps) 

mmm rr тортны, RR RU ВАС 8 А ЬЬ НЫ Ч ое 2 AA CORR EEE · 
1.0247 -8.609 1.849 
1.0254 -10.380 1.816 
1.0260 -11.528 1.783 
1.0267 -11.982 1,751 
1.0273 -11.716 1.718 
1.0280 -10.745 1.684 
1.0286 -9.129 ai 
1.0293 -6.964 1.617 
1.0299 -4.380 1.584 
1.0306 -1.534 1.550 
1035183 1.404 1.516 
1.0319 4.258 1.482 
1.0326 6.855 1.447 
1.0332 9.040 1.413 
1.0339 10.682 1.378 
1.0345 11.681 1.344 
150352 11.977 1.309 
1.0358 11.554 1.274 
1.0365 10.435 1.239 
1.0371 8.689 1.204 
1.0878 6.421 1.168 
1.0384 3.766 1.133 
1.0391 0.885 1.097 
1.0397 -2.050 1.062 
1.0404 -4.862 1.026 
1.0110 -7.382 0.990 
1.0417 -9.458 0.954 
1.0423 -10.966 0.918 
1.0430 -11.815 0.882 
1.0436 -11.953 0.846 
1.0433 =] 1373 0.810 
1.0449 -10.109 0 773 
1.0456 -8.237 0.737 
1.0462 -5.869 0.700 
1.0469 -3.147 0.664 
1.0475 -0.236 0.627 
1.0482 2.691 0.591 
1.0488 5.456 0.554 
1.0495 7.898 0.517 
1.0501 9.859 0.480 
1.0508 11.232 0.443 
1.0514 11.930 0.406 
1.0521 11.910 0.369 
1.0527 IU PTS 02932 
1.0534 9.768 0.295 
1.0540 7103 0.258 


Table B-2. (cont.) Sample simulation output file. 
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ТТМЕ (5ес) 


. 0547 
·. 0553 
.0560 
.0566 
.0573 
. 0579 
.0586 
‚0592 
‚0599 
‚0605 
‚0612 
‚0618 
‚0625 
- 0632 
‚0638 
‚0645 
‚0651 
‚0658 
‚0664 
‚0671 
‚0677 
‚0684 
‚0690 
‚0697 
‚0703 
‚0710 
‚0716 
‚0723 
. 0729 
.0736 
‚0742 
‚0749 
20755 
.0762 
.0768 
‚0775 
‚0781 
‚0788 
‚0794 
1.0801 


ka іі ыі іі кі кӛ мі іі кі і-4 ә sas as eA fateh fh fm Іні ә ма һә ма а fh fee fee а а Бә pe 


Table B-2. 


STATOR PHASE A ROTOR LOOP 1 
CURRENT (amps) CURRENT (amps) 
5.310 0.221 

2 526 0.184 
-0.410 0.147 
-3.324 0.110 
-6.038 0.072 
-8.392 0.035 
-10.242 -0.002 
-11.477 -0.039 
-12.024 -0.076 
-11.848 -0.113 
-10.960 -0.151 
-9.413 -0.188 
-7.300 -0.225 
-4.746 -0.262 
-1.906 -0.299 
1.050 -0.336 
3.945 -0.373 
6.604 -0.410 
8.869 -0.447 
10.602 -0.484 
11.700 -0.521 
12.095 -0.558 
11-765 -0.594 
10.727 -0.631 
9.045 -0.668 
6.819 -0.704 
4.181 -0.741 
1.291 -0.777 
-1.679 -0.813 
-4.551 -0.850 
-7.151 -0.886 
-9.323 -0.922 
-10.938 -0.958 
-11.897 -0.994 
-12.143 -1.030 
-11.660 -1.066 
-10.477 -1.101 
-8.665 -1.137 
-6.331 -1.172 
-3.616 -1.207 


(cont.) Sample simulation output file. 
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TIME (sec) STATOR PHASE AÀ ROTOR LOOP 1 


CURRENT (amps) CURRENT (amps) 
1.0807 -0.681 -1.243 
1.0814 2.297 -1.278 
1.0820 5.139 -1.313 
1.0827 7.675 -1.347 
1.0833 9752 -1.382 
1.0840 11.246 -1.417 
1.0846 12.067 -1.451 
1.0853 128511615 -1.485 
1.0859 181825 3/3 -1.519 
1.0866 10.209 -1.553 
1.0872 8.272 -1.587 
1.0878 5.838 -1.621 
1.0885 3.051 -1.654 
1.0892 0.080 -1.688 
1.0898 -2.899 -1.721 
1.0905 -5.706 -1.754 
1.0911 -8.173 -1.787 
1.0918 -10.152 -1.820 
1.0924 -11.524 -1.852 
1.0931 -12.207 -1.884 
1.0938 -12.158 -1.917 
1.0944 -11.381 -1.949 


Table B-2. (cont.) Sample simulation output file. 
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STATOR PHASE CURRENT 


SIDE 


SIMULATION 


04 1728 rpm 


ROTOR BAR 


ONE BROKEN 


балку 





UL ІМЯЧЯПӘ 


1.133 1.200 
TIME in SECS 





ОНИ МЕ ТОМУ ROTOR LOOP CURRENT 
EXER s DEIN B rpm 
ONE BROKEN ROTOR BAR 


CURRENT in AMPS 





1.000 ISS 1.267 1.400 us 1.667 
TENE in SEES 


Figure B-2. Rotor loop 1 current vs. time for sample 
simulation. 
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B-2 PRO-MATLAB EIGENVALUE ROUTINE 


PRO-MATLAB is a matrix computation software package 
developed by The MathWorks, Inc. of Sherborn, MA. This 
PRO-MATLAB script file calculates the eigenvalues of the 
[RR]/[L] matrix. The eigenvalues are computed for rotor 
slipz1.0 (i.e., the rotor is stationary). 


a. PROGRAM LISTING 
program EIGENVALUES 


filename: root.m 
author: Mark S. Welsh 
references: PRO-MATLAB user's guide 


PRO-MATLAB routine to calculate eigenvalues for induction 
motor system of equations ([r]/[1]). 


variable list 
sf-rotor slip 
t=time 
nrb=no. of rotor bars 
rleq,11eq,112eq,12eq,r2eq-eq. circuit model values 
pp=machine pole-pairs 
wzinput electrical frequency 
wm=mechanical speed 
k=constant 
lsl=stator phase leakage inductance 
rs-stator phase resistance 
ls=stator phase self inductance 
m=mutual inductance coefficient for stator-rotor 
lrl=rotor bar leakage inductance 
rrzrotor bar resistance 
lr=rotor loop self inductance 
1,ll=inductance matrix 
r,rlzresistance matrix 
meigzeigenvalues 
tl=smallest time constant 
t2=zlongest time constant 


26 36 96 96 96 92 92 56 56 56 92 56 C IC JE ae JL Be эё JE BE IC IE FE Be ae at at ae IE ae Ie 


input parameters 


sf=0; 

t=0; 

nrb=45; 
гјед=. 859; 
lleq=.0046; 
112eq=.0704; 
12eq=.0046; 
r2eq-.5612; 
pp=2; 

EK. 


X calculate electrical and mechanical parameters 
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wm=wt (l-sf)/pp; 
k=pi*pp/nrb; 


Шеш 


lleq; 


rs=rleg; 
1s=2*112eq/3; 
mz8t*sin(k)*112eq/(3*pi); 


1г1= 


(12еа-(К^2/(вї1їп(К)^2)-1)%112ед)%4%пгЬ/(3%р1^2); 


гг=4%пгЪ/(3®р1^2)*г2еа; 
lrz16t(nrb-1)*pp^2*112eq/(3tnrb^2); 

* calculate elements of r and 1 matrix 
fi EE MESE 


end 
for 


end 


1(4, 
r(4, 


for 


end 
% reduce 
for 


HO Э 
r(i,j)=:0; 
1(i,j)=-1s/2; 
end 
ri тв 
lti,i)zls#+1s]; 
for 3J=4:nrb+3; 
l(i,j)zm*sin(pptwm*t*«(2t(j-3)-1)*pi*pp/nrb-(i-1)* 
2*pi/3): 
r(i,j)-pp*wmtmtcos(pp*wmtt*(2*(j-3)-1)*pitpp/nrb- 
(1-1)%2%рі/3); 


епа 
1=4:пгЪ+3; 
for DS 
ТА = O 
1(1, ј)=-1г/(пгб-1); 
if ј==1+1: 
l(i,j)--lr/(nrb-1)-1rl; 
r(i,;))--rr; 
end 
end 


Tel (Mle i =o trr: 
(а) = Ота и 


nrb+3)=-1r/(nrb-1)-1rl; 
nrb+3)=-rr; 
1=1:пг6+3; 
for j-i:nrb*3; 
TIGUPEIBIE r ( 0) 
ПО Е) 
епа 


system for sum of rotor loop currents=0 
i-1:nrb*2; 
for jz1:nrb*2; 
if 3€4; 
ia.) =D) 
else; 
r1l(1i,)j)=r(i,))-r(i,nrb+3); 
11(1,23):1(13,2)-1(1,nrb*3); 
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епа 
епа 
епа 
calculate eigenvalues 
meigz=eig(-ll\rl); 
determine smallest time constant 
tl=1/max(abs(meig)); 
determine largest time constant 
t2=1/min(abs(meig)); 
save values 
save eigs.mat; 
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b. EIGENVALUES FOR 3-BAR ROTOR 


The following values are the eigenvalues calculated for 
the 3-bar rotor discussed in Chapter 3. 


-186.740 
-154.590 
-154.590 
-4.662 
-4.662 


Table B-3. 3-bar rotor eigenvalues. 
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c. EIGENVALUES FOR 45-BAR ROTOR 


The following values are the eigenvalues calculated for 


the 45-bar rotor (1.е., experimental motor) discussed in 
Chapter 3. 
-186.740 -80.344 
-154.590 -80.344 
-4.662 -84.698 
-154.590 -88.364 
-4.662 -88.364 
-1.961 -93.992 
-1.961 -91.436 
-15.636 -91.436 
-15.636 -93.992 
-25.273 -96.101 
-25.273 -96.101 
-35.354 -97.817 
-35.354 -97.817 
-45.124 -99.185 
-45.124 -99.185 
-54.129 -101.010 
-54.129 -101.510 
-62.161 -101.010 
-62.161 -101.510 
-69.171 -100.240 
-69.171 -100.240 
-75.203 -101.760 
-75.203 -101.760 
-84.698 


Table B-4. 45-bar rotor eigenvalues. 


198 





B-3 PRO-MATLAB FFT ROUTINE 


This PRO-MATLAB script file computes the fast Fourier 
transform of the stator phase and rotor loop current 
outputs of the FORTRAN simulation program. The magnitude of 
the FFT values are converted to dB and written to a file for 
future plotting if desired. 


a. PROGRAM LISTING 
program FFT 


filename: fast.m 

author: Mark S. Welsh 

references: PRO-MATLAB user's guide 
references [18] ,[19] 


PRO-MATLAB procedure to calculate the FFT of the stator 
phase and rotor loop current outputs from the FORTRAN 
simulation program. The input data is sampled at a 
frequency of 153.6 Hz (or 512 data points are taken). 

A second-order Hanning window is used to aid in 
resolution of harmonic frequency components present. 
The power spectral density function is computed and 
converted to a dB scale. | 


variable list 
a-input signals array (t,ia,irl) 
t-time 
ia-stator current 
ir=rotor loop 1 current 
s=FFT of ia 
r=FFT of ir 
ps=power spectral density of ia 
pr=power spectral density of ir 
w=hanning weight factors 


IL IL IL JL JL FL JE JL JE JL FL JL JL JE IL IL FL зе JL IL IL IL FL IL IL al ae ae 


sample data at 153.6 hz 
t=a(1.10:5120,1); 
pasa (1: 105120, 2) ; 
іт-а(1:10:5120,3); 
% calculate window function 
wz.55(0ones(1:512)-cos(2*pi*(0:511)/512)); 
X weight inputs using window function 
ijiawzia.tw'; 
ШЕ. Аы: 
* determine ffts 
s=FFT(iaw) ; 
r=FFT(1rw); 
% determine spectral densities 
ps=s.%conj(s); 
przr.*conj(r); 
* convert to dB 
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ps:10*log10(ps); 
pr=10*log10(pr); 

% determine frequencies 
ISS A OOO 502; 

% save values and put in plot format 
bell ps( 12256) ); 
esp wbr(i1:256)]; 
save ffts.mat; 
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b. SAMPLE FFT OUTPUT FI 


The following table is the output from the PRO-MATLAB 
FFT routine. The values given are the power spectral 
densities for the FORTRAN simulation program output currents 
of section 1.1.2 above. Plots of the frequency spectrum for 
stator phase a and rotor loop 1 currents are included. 





FREQUENCY STATOR PHASE A ROTOR LOOP 1 
(Hz) CURRENT PSD CURRENT PSD 
(dB) (dB) 

0.0 -120.490 -48.451 
0.3 -117.140 -53.454 
0.6 -121.350 -68.221 
0.9 -117.370 -63.254 
И 2 -117.610 -56.369 
155 -119.050 -48.357 
1.8 -115.980 -36.333 
2 o Jl -120.710 47.687 
2.4 -118.420 53.710 
2o -117.390 41.692 
3.0 -120.740 -36.337 
3.3 -120.190 -48.388 
3.6 -117.220 -56.358 
229 -118.840 -62.401 
4.2 -117.980 -67.268 
4.5 -120.120 -71.393 
4.8 -117.740 -74.918 
В -120.040 -78.041 
5.4 -118.100 -80.848 
Sn -118.700 -83.388 
6.0 -117.840 -85.617 
BS -118.500 -87.782 
6.6 -116.340 -89.849 
6.9 -119.570 -91.454 
7.2 -117.290 -93.317 
1.9 -118.280 -94.946 
7.8 -115.270 -96.140 
8.1 -119.600 -98.164 
8.4 -116.090 -98.824 
8.7 -116.310 -100.810 
9.0 -119.770 -101.670 
3,3 -113.980 -102.490 
326 -118.430 -104.630 
229 -116.180 -104.240 
10.2 -117.560 -106.670 
10.5 -115.430 -106.870 


Table B-5. FFT output file. 
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FREQUENCY STATOR PHASE A ROTOR LOOP 1 


(Hz) CURRENT PSD CURRENT PSD 
(dB) (dB) 

AAA ——— P  —T-.M. —=— = — í —— r — ee LLLA — 
10.8 2157520 -107.940 
ШЕН! -116.440 -109.460 
11.4 -114.510 -108.920 
11.7 -118.390 -111.340 
12.0 -114.800 -112.150 
12,3 -113.590 -110.150 
12.6 -116.000 -116.330 
12,3 -114.220 -111.390 
13772 -114.280 -115.520 
IE -116.640 -115.190 
13,8 -113.160 -115.720 
14.1 -113.980 -115.860 
14.4 -114.900 -119.660 
14.7 -114.180 -118.010 
15.0 -113.710 -116.340 
123 -113.610 -118.070 
15.6 -113.860 -125.450 
oro -113.600 -119.010 
16.2 -111.620 -117.840 
16.5 -113.520 -124.520 
16.8 -112.260 -120.970 
wW. l -112.840 -121.440 
17.4 -112.180 -122.970 
NT T -112.190 . -124.910 
18.0 -112.020 -121.540 
18.3 -111.010 -125.780 
18.6 -112.460 -120.890 
18.9 -110.760 -120.280 
23,2 -112.040 -115.820 
19.5 -110.670 -124.750 
19,8 -109.900 -121.160 
20.1 -111.180 -125.260 
20.4 -111.230 -124.470 
20.7 -108.920 -142.480 
21.0 -110.990 -122.410 
21.3 -99.670 -131.490 
21.6 -98.796 -123.640 
21-9 -99.621 -127.000 
22-2 -109.030 -129.970 
22.5 -109.380 -128.540 
22.8 -108.480 -125.770 
2351 -108.520 -126.560 
23.4 -107.810 -135.920 
23.1 -108.910 -125.290 


Table B-5. (cont.) Sample FFT output file. 
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ЕКЕОЏЕМСУ STATOR PHASE A ROTOR LOOP 1 


(Hz) CURRENT PSD CURRENT PSD 

(dB) (dB) 
24.0 -107.210 -129.800 
24.3 -107.170 -127.430 
24.6 -108.180 -132.510 
24.9 -107.020 -130.390 
25,2 -106.790 -133.600 
49.5 -106.530 -129.600 
25.8 -107.180 -124.460 
26.1 -105.870 -126.050 
26.4 -105.780 -125.480 
Ст -105.290 -131.720 
24 .0 -106.050 -132.720 
27.3 -105.470 -140.190 
27.6 -104.450 -122.270 
27.9 -105.650 -129.900 
28.2 -104.010 -135.800 
28.5 -104.180 -125.730 
28.8 -104.100 -134.130 
29.1 -103.810 -134.500 
29.4 -103.620 -124.620 
29.7 -102.930 -125.320 
30.0 -103.550 -126.870 
30,3 -102.190 -125.760 
30.6 -102.490 -119.520 
30.9 -102.120 | -123.740 
31.2 -102.120 -129.180 
21-5 -101.730 -124.990 
31.8 -101.060 -132.810 
22.1 -101.390 -135.670 
32.4 -100.340 -129.790 
2.7 -100.610 i -126.290 
33.0 -100.460 -132.600 
99.3 -99.795 -130.500 
33.6 -99.354 -126.870 
99.9 -99.584 -137.510 
34.2 -98.968 -136.360 
34.5 -98.627 -144.390 
34.8 -98.221 -131.530 
95.1 -98.288 -145.710 
35.4 -97.658 -126.560 
SoG? -97.524 -76.744 
36.0 -96.991 -70.852 
36.3 -96.985 -76.980 
36.6 -96.150 -128.330 
36.9 -96.247 -128.480 
97.2 -95.819 -124.140 


Table B-5. (cont.) Sample FFT output file. 
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FREQUENCY STATOR PHASE A ROTOR LOOP 1 


(Hz) CURRENT PSD CURRENT PSD 

(dB) (dB) 
37.5 УВО -129.900 
37.8 -95.145 -148.670 
38.1 -94.592 -128.070 
38.4 -94.500 -131.560 
Sorat -94.062 -129.210 
39.0 -93.424 -130.500 
39.3 -93.278 -123.870 
39.6 -92.801 -133.130 
39,9 -92.563 -127.250 
40.2 -92.018 -136.670 
40.5 -91.743 -110.170 
40.8 -91.323 -104.290 
41.1 -90.756 -110.790 
41.4 -90.485 -126.750 
41.7 -89.964 -133.240 
42.0 -89.702 -128.000 
42.3 -89.106 -129.450 
42.6 -88.737 -129.520 
42.9 -88.230 -133.660 
43.2 -87.867 -139.180 
43.5 -87.242 -128.160 
43.8 -86.899 -144.760 
44.1 -86.351 -127.010 
44.4 -85.971 -127.780 
44.7 -85.277 -129.190 
45.0 -84.879 -135.890 
45.3 - 84.398 -140.460 
45.6 -83.763 -128.750 
45.9 -83.270 ` -128.910 
46.2 -82.767 -132.750 
46.5 -82.112 -132.010 
46.8 -81.559 -131.580 
47.1 -81.015 -131.070 
47.4 -80.360 -138.210 
47.7 -79.759 -125.470 
48.0 -79.152 -129.590 
48.3 -78.467 -137.420 
48.6 -77.825 -130.090 
48.9 -77.168 -135.140 
49.2 -76.467 -132.980 
49.5 -75.738 -131.590 
49.8 -75.034 -136.100 
50.1 -74.285 -148.390 
50.4 -73.528 -136.750 
50.7 - 72.760 -131.330 


Table B-5. (cont.) Sample FFT output file. 
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FREQUENCY STATOR PHASE A ROTOR LOOP 1 


(Hz) CURRENT PSD CURRENT PSD 

(dB) (dB) 
51.0 -71.924 -141.820 
51.3 -71.151 -136.470 
51.6 -70.303 -157.810 
51.9 -69.445 -133.010 
52.2 -68.565 -129.240 
52.5 -67.630 -132.530 
52.8 -66.607 -127.240 
53.1 -65.291 -134.450 
53.4 -63.277 -131.820 
53.7 -59.719 -136.240 
54.0 -53.913 -134.880 
54.3 -45.531 -130.020 
54.6 -33.052 -128.300 
54.9 23.979 -131.300 
55.2 30.056 -134.000 
55.5 24.092 -135.190 
55.8 -32.168 -123.950 
56.1 -41.944 -124.990 
56.4 -44.745 -130.200 
56.7 -44.090 -121.450 
57.0 -42.096 |^ — -123.830 
57.3 -39.464 -126.810 
57.6 -36.775 -123.950 
57.9 -33.185 -125.870 
58.2 -29.069 -132.570 
58.5 -24.232 -138.080 
58.8 -18.225 -134.530 
59.1 -10.281 -134.520 
59.4 1.728 -133.890 
59.7 57.809 -129.160 
60.0 63.892 -135.130 
60.3 57.932 -137.680 
60.6 1.878 -129.910 
60.9 -10.194 -131.000 
@1.2 -18.165 -129.220 
61.5 -24.193 -127.560 
61.8 -29.059 -126.040 
Є2 1 -33.145 -123.560 
62.4 -36.669 -134.160 
62.7 -39.769 -132.280 
63.0 -42.536 -136:570 
63.3 -45.037 -128.000 
63.6 -47.316 -124.010 
63.9 -49.412 -127.040 
64.2 -51.350 -134.290 


Table B-5. (cont.) Sample FFT output file. 
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FREQUENCY STATOR PHASE A ROTOR LOOP 1 


(Hz) CURRENT PSD CURRENT PSD 

(dB) (dB) 
64.5 -53.152 -130.420 
64.8 -54,841 -130.330 
65.1 -56.422 -133.160 
65.4 -57.911 -145.320 
69.7 -59.326 -135.670 
66.0 -60.656 -125.890 
66.3 -61.930 -127.300 
66.6 -63.138 -139.560 
66.9 -64.300 -128.220 
67.2 -65.392 -132.960 
СЭ -66.459 -125.600 
67.8 -67.466 -127.090 
68.1 -68.455 -128.730 
68.4 -69.369 -143.340 
68.7 -70.296 -135.280 
69.0 -71.153 -128.240 
69.3 -71.998 -128.200 
69.6 =Y 2o YNZ -133.360 
69.9 -73.598 -127.150 
70.2 ~74.317 =121.150 
noc -75.079 -123.520 
70.8 =75.729 -128.960 
ШІСІ -76.456 -137.850 
71.4 -77.050 -149.120 
el of -77.736 -134.230 
72.0 -78.304 -132.220 
72.3 -78.885 -135.290 
72.6 -79.435 -139.460 
299 -79.952 -126.330 
23.2 -80.441 -127.550 
73.95 -80.928 -131.900 
7328 -81.307 -129.540 
74.1 -81.752 -136.740 
74.4 -82.122 -129.350 
74.7 -82.474 -132.450 
75.0 -82.672 -131.340 


Table B-5. (cont.) Sample FFT output file. 
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ПО О АУЕ ОДО УВ ОДЕ ЕРИ SPESTRUN 
SLIP=.04 1728 rpm 
ONE BROKEN BAR 


in dB 


SPECTRAL POWER DENSITY 





50 52 54 536 58 60 62 64 66 68 70 
Blogos C Y im HZ 


Figure B-3. Stator phase a frequency spectrum for sample 
simulation. 
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ROTORZEOOEZEFREOUENEY SPECTRUM 
SLIP=.04 1728 rpm 
ONE BROKEN BAR 


ШЕ «ІН 


POWER DENSITY 


юш ЕЛҮ 





DUCIT No а о 267,187 20 
БББ im EZ 


Figure B-4. Rotor loop 1 frequency spectrum for sample 
simulation. 
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APPENDIX C 


EXPERIMENTAL DATA 


The following sections contain tables of the measured 
data (the fundamental and (1-2s)f components of the stator 
phase current) for the four rotors tested. This is the data 
recorded using the HP-3561A dynamic spectrum analyzer as 
discussed in section 4.2. Included in each table is the 
average value and standard deviations for each data set. In 
addition, the average stator phase voltage is included in 
the tables for the fundamental component of the stator phase 
current. An analysis of the data presented below is included 
in Chapter 4. 
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C-1 EXPERIMENTAL DATA ROTOR #1 














КОТОК #1 
I(1-2s)f COMPONENT OF STATOR CURRENT 
AMPS rms 


0.0014 0.0026 0.0027 
2 0.0010 0.0017 0.0019 0.0034 






[0.0012 | 0.0028 | 0.0032 | 0.0027 | 
s [0.0008 | 0.0018 | 0.0020 | 0.0040 | 
[0.0011 | 0.0017 | 0.0036 | 0.0017 _| 
7 [0.0008 [0.0016 | 0.0019 | 0.0001 | 
[0.0007 | 0.0012 | 0.0018 | 0.0080 _| 
[0.0012 | 0.0013 | 0.0019 | 0.002 | 
[0.0008 | 0.0011 | 0.0028 | 0.0037 | 
6.0010 | 0.0018 | 0.0024] 0.0018 | 
[0.0007 | о.00 | 0.0019 | 0.0040 | 









[ 6.001 | 


0.0012 0.0015 0.0024 0.0020 


” EA 
кі 


10 
12 
13 
14 
15 
16 
167 






MEAN 0.0011 | 0.0015 | 0.0025 | 0.0032 
STD DEV. 0.0003 0.0005 0.0005 0.0009 | 


Table. C-1. ROTOR #1 I(1-2s)f raw data. 
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ROTOR #1 
FUNDAMENTAL COMPONENT OF STATOR CURRENT 
AMPS rms 


ТЕТЕ owe Ж !erence] 097 | 
eae Ae a 
waris 2.02 | 

8.98 


ји a ueri. 


AVG INPUT 
123.3 123.1 TNT 122.5 | 


VOLTAGE 
Table C-2. ROTOR #1 1(60 Hz) raw data. 




























(V rms) 





C-2 EXPERIMENTAL DATA ROTOR #2 (ONE END OPEN) 






ROTOR #2 
I(1-2s)f COMPONENT OF STATOR CURRENT 
AMPS rms 


Tiesto, [ ser T 3:9.» [ 3:939  [ 2:923 ] 
[1] 0.0020] 0.0018 | 9:00: | 0.0012 | 
[— 3 | 901: | oon [ 903 | 009 ] 
[0.0011 | 0.0014 | 0.0020 | 9.001: | 
Ts] 0.0008 | 0.0028 | 0.0026 | 0.0029 | 
0.0013 | 0.0032 | 0.0009 “| 9.0036 | 
[7 [0.0015 | ooon | 0.0028 | 0.0020 | 
[I [0.0015 | 0.0013 | 0.0020 | 0.0082 | 
[э | 0.0010 | 0.0019 | 0.0022 | 0.0015 _| 
[10 [0.0010 | 0.0018 | 0.0029 _| 0.0017 | 
E 


Table C-3. ROTOR #2 I(1-2s)f raw data. 



























р 
[ 
| 
н 


| 
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ROTOR #2 
FUNDAMENTAL COMPONENT OF STATOR CURRENT 
AMPS rms 


шоп |] 
аи 
as 


MEAN 
STD DEV. 


AVG INPUT 
VOLTAGE 
(V rms) 





Table C-4. ROTOR #2 1(60 Hz) raw data. 
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C-3 EXPERIMENTAL DATA ROTOR #2.1 (BOTH ENDS OPEN) 


ROTOR #2.1 
I(1-2s)f COMPONENT OF STATOR CURRENT 
AMPS rms 


тыт [ший [50 [оз [тоот 
[1] 0,008 | 0.0083] 0.0101] 0.0169] 
[a so | 0.0070 | 0.0121 | 9.0184 | 
ps [e-0056 [0.0067 | 0.0120 | 0.0175 
[= [0.0077 | 0.0075 | 0.0128 | 0.0198 | 
т [| sos | 0.0078 | 0.0181 | 0.0185 | 
[s 0.0088 | 0.0087 | 0.0147 | 0.0188 | 
| | een | оов [сёз [ое | 
NN Баста голе u |o mi 
Tess | woe 1 0-08 | ова | 
| 0.0145 j 0145 


Table C-5. ROTOR $2.1 I(1-2s)f raw data. 
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ROTOR #2.1 
FUNDAMENTAL COMPONENT OF STATOR CURRENT 
AMPS rms 


Lr e — 


Sesh es) oe) m sme | 
[s [oes | ыы | та еле | 


Fen | ам | вез Тын” | њег | 


AVG INPUT 
VOLTAGE 123.1 122.9 122.5 122.2 
(V rms) 





Table C-6. ROTOR #2.1 I(60 Hz) raw data. 
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C-4 EXPERIMENTAL DATA ROTOR ¿3 









ROTOR #3 
I(1-2s)f COMPONENT OF STATOR CURRENT 
AMPS rms 


rest no] 820.01] 20.02] 820.08 [оо 
[1 | 0.0008 | omo [ 0.0011 | 0.0019] 
[2 [0.0007 | 0.0008 | 0.0014 | 9.0013 — 















0.0005 0.0006 | 0.0012 | 0.0016 
5 0.0009 0.0006 0.0011 0.0016 





[s 71 0.0005 | 0.0008 | 0.0018 | 0.0012 | 
[s | 0.0006 | 0.0008 [0.0012 | 0.0011 _| 
[100-0005 | 0.0008 | 0.0011 | 0.0018 


MEAN 

























^. 0.0007 0.0006 ` 0.0011 | 0.0013 | 
STD DEV. 0.0003 0.0002 0.0003 0.0002 | 


Table. C-7. ROTOR £$3 I(1-2s)f raw data. 
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КОТОК #3 
FUNDAMENTAL COMPONENT OF STATOR CURRENT 
AMPS rms 


[test no. | 20.01] 70.02 | s0.09 Гон 
Peso seo tao | eer | 
ee 
_ ma sa вука | теле есе | 
[a o; p $5 p ГГ | 
Im n e su = = 
Ee e ens 
MEL a a 


AVG INPUT 
1225 7 122.4 122.1 121.8 


VOLTAGE 
Table C-8. ROTOR #3 1(60 Hz) raw data. 


























(V rms) 
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APPENDIX D 


INDUCTION MOTOR STATOR PHASE CURRENT HARMONICS 


D-1 Introduction 


In addition to the fundamental and (1-2s)f frequency 
components of the stator phase current, components at 
frequencies (1-4s)f, (1+2s)f, and (1%4а)Ғ меге а1во 
detected for each rotor tested as shown in Chapter 4. In 
order to examine the source of these additional frequency 
components, an expression for the frequency components 
which can exist in the stator phase currents due to both 
time- and space-harmonic components of the airgap flux is 
derived. 


D-2 Derivation 


Assume initially that the airgap flux density due to 
currents flowing in the stator phase windings is given by 
the following expression 


B. = Е Д {B,,sin(mwt+npé@) 


m-oddt-odd 


+ B,,sin(mwt-np0)) (D-1) 


where 


B,,.B,,* Fourier coefficents 
m*m "stator time-harmonic 


n=n"stator space-harmonic 


Using the transformation relation given by equation 
2-5, in the rotor reference frame equation D-1 becomes 


pru > » (B,,sin(mwot*np80'*n(l-s)ot) 


m-eddn-odd 


-B,.sin(mwt-npd°-n(l-s)wt)) (D-2) 


218 





From equation D-2, it can be seen that airgap flux will 


couple with the rotor loops and induce currents of 
frequencies 


(m= n(1-s))O for m,n odd (D-3) 


These induced rotor currents also produce an airgap 
flux which can be expressed in the following form 


В,= ) ) ) {Byassin((m+n(1-s))w t+ qpe’) 


a-odd m -oddn-edd 
+ В сь 510 ((т + п(1-5))01-4рӨ”) 
= Ела и (а = п(-—=))о ар) 


* Bm sin((m-n(1l-s))wt-qp9’)} (D-4) 


where 


B Sp onp Ba." bourrer coefficients 


q= q" rotor space-harmonic 


In the stator reference frame this becomes 


B.* ) 2 2 «Ви sin((men(I-8)- q(1-s))wt+qp0) 
q-04¢ m-0ddn-04d 
+ Bae, sin((m~+ n(i-s)* q(1-s))a@t-qpée) 
+B. sin((m-n(I-s)-q(1I-s))eot+qp0) 
+ Bynosin((m-n(l-s)+ q(!l-s))wt-qp0)) (D-5) 
This component of the airgap flux will couple with the 


stator phase windings and induce currents with the 
following frequencies 
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(m*n(1-s)-q(1-s))w for m.n.q odd (ето 


(m-+-n(1-s)+ q(1-s))O for m,n,q odd OD 70 
(m-ní(l-s)-q(1-s))o for m.n.q odd (DEBI) 
(т-п(1-5) + 4(1-5)) о for m,n,q odd (D-9) 


Thus for each combination of m, n, and q there can be 
as many as four distinct frequencies of currents induced in 
a stator phase winding. Table D-1 shows the possible stator 
phase current harmonics for m,n,q=1,3,5. This table shows 
that the (1-4s)f and the (1+2s)f components are produced by 
the third time-harmonic and that the (1+4s)f component is 
produced by the fifth time-harmonic component of the stator 
phase current. 


For the experimental motor with ROTOR #3 installed and 
operating at a slip of 0.0011, the stator phase current 
frequency spectrum measured using the HP-3561A signal 
analyzer is shown in Figure D-1. This spectrum shows that 
the magnitudes of the third and fifth time-harmonics are 
approximately the same and that they are larger than all 
other harmonic components except the fundamental component 
of the stator phase current. This would indicate that for 
the motors tested, the (1-4s)f, (1+2s)f, .and (1+4s)f 
components can be attributed to the third and fifth 
time-harmonics of the stator phase current. 
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time/space STATOR PHASE CURRENT HARMONICS 
component normalized to the fundamental frequency 
(eqn D-6) (eqn D-7) (eqn D-8) 
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Table D-1. Predicted stator current harmonics. 
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Figure D-1. Stator current time-harmonics. 
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